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(54) TiUe: MULTICOLOR OPTICAL PYROMETER 
(57) Abstract 

A method of using multiple single-color pyrometers 
to correct for die effects of reflected radiation when optically 
measuring the temperature of a relatively cool surface in 
hotter surroundings is disclosed. The method includes 
obtaining pyrometric readings from the surface at two 
or more wavelengths and thereafter utilizing tiie readings 
together with the emissivity of the surface to calculate the 
surface temperature. In one version of the method the 
emissivity of the surface being measured is deterrnined 
by a method which includes heating the surface in an 
insulated environment and thereafter obtaining readings 
pyrometrically from the surface once the insulation of the 
surface is removed. The readings are then utilized to 
calculate the emissivity by extr^lating the readings to a 
time zero. In those versions' of the method which utilze 
pyrometric readings taken over three or more wavelengths, 
emissivity may be calculated through the use of relevant 
mathematical expressions. 
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MULTICOLOR OPTICAL PYROMETER 



BACKGROUND OF THE INVENTION 
5 Technical Field : The invention relates lo optical measurement. The 

invention is specifically directed to the measurement of temperatures using infrared 
technology. 

Statement of tlie Art : Tliere ^are many applications of radiation thermometers 
throughout the industrial world. The first modem attempt at radiation thermometry 

10 dates as far back as 1828 (Dike et aL, 1966). The first comprehensive treatment of 
the theory of radiation thermometry and its practice was published in 1960 by 
Thomas R. Harrison (1960). Since then, many of the most important advances in 
this rleld have been described and chronicled in volume 4 and 5 of the well-known 
proceedings, Temperature, its Measurement and Control in Science and Industry, a 

15 series of symposia beginning in 1929. More recently, DeWitt and Nuuer (1988) 
have provided a comprehensive coverage of the principles and applications of 
radiometric methods for measuring surface temperature. 

The following U.S. Patents have been identified as identifying representative 
effons in this area: 3,619,059; 3,454,769; 4,020,695; 4,818,102; 5,029,117; 

20 5,125,739; 5,226,732; 5,231,595; 5,318,362; 5,326,173. 

Ono (1988) and Nutter (1972) cover a number of previous investigations to 
cope with the extraneous radiation caused by hotter surroundings and with uncertain 
emissivities. 

Atkinson and Strange (1976) developed a technique using a fiber pyrometer 
25 probe with two identical silicon diode detectors to measure the temperature of 

turbine blades in the presence of reflected energy originating from the flame. The 
spectral bands they chose were 0.4 -1.2 |tm (as received or unfiltered) and 0.4-0,85 
^^m (filtered). Note that the latter range overlaps with the former. They described 
their system with the following equations 
30 = llsX^)P[T,)dX + cj;5,(A)P(rJdA 

(A) 

^.,o.u=r,S,{X)P{T,)dX^C^\''^S,{X)P{T^)dX 
where E, is the reading of the pyrometer, and T, are the blade temperature and 
the reflected source temperature; Sj(X) is the sensitivity of pyrometer i at 
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wavelength X; P(T) is the monochromatic emissive power of a blackbody ai 
temperature T and C, is a complex, constant term that includes the effects of the 
emissivity of the surface, the emissivity of the gas, and geometiacai considerations 
of reflected radiation. 

5 The investigators assumed Si(X) could be pre-determined from calibration 

studies and from knowing the spectral sensitivity of the detectors. They further 
assumed T, was measured or was estimated from the general characteristics of 
flames in turbine engines. Then only two variables, namely Tb and C^, were left as 
unknowns in the two expressions of equation (A). A trial and error method 
10 involving graphical analysis was used in an attempt at obtaining the blade 
temperature. 

Some of the assumptions in this study are questionable. The authors treated 
the turbine blade as a blackbody on the basis of the nearly isothermal conditions 
within the engine (therefore no emissivity factors were needed in the first terms of 
15 the equation (A) relating to blade temperature, TJ. An extremely hot flame 

2148'='C (3900 **F) was assumed as the reflected energy source and the temperature 
difference between the leading edge of the blade and body of the blade was taken as 
220'*C (400''F). The reliability of using the same Si(X) for both T^, and T, is also 
doubtful. 

20 DeWitt (1986) reviewed and commented on the similar work by Imgram and 

McCanless (1983) which aimed at determining the surface temperature of tubes in a 
crude oil refinery heater. By sighting through a peephole, the radiation 
thermometer viewed a portion of the tube, and the problem was to relate the 
pyrometer signals to the true tube wall surface temperature. The mathematical 

25 system was given as 

£«(7'«) = f.:£.:..(T,) + F_(l-£,,)E,,,(r^,) + (l-F_)(l-£,,)£,^,(r,,) 

where is the view factor between the furnace wall and the target surface, T^j is 
30 the wall temperature of the funiace, and Tn is the temperature of the flame. T^j and 
Tn were determined optically. The emissivities, e^. and e>^, were known from other 
measurements. 
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Waveiengths of 0.9 ^tm and 3.9 were used and the emissivity of the 
surface of the tube was 0.9 for both and X.. Table 1 shows the results obtained 
(Dewitl, 1986). Tlie T^'s and T^'s were measured optically at 0.9 and 3.9 fim. 
Tube surface temperatures, Tx;, which included the reflected components, were also 
) measured in the same manner at two locations (one was closer to the flame and the 
other is closer to the wall). The true surface temperature of the tubes, T„ and the 
view factor, F^.„ were calculated from equations (B). The differences between the 
optically measured temperatures and the calculated temperatures, T^i-T,, in different 
locations and at different wavelengths, are also tabulated. 

Table 1 Measurements in a direct-fired, crude oil refuiery heater (Dewitt, 1986). 

T^i-^C Tx, . Ts, °C 



15 



30 



0.9 ^m 3.9 Tw.t 0-9 l^m 3.9 

name Temp. Tf 1521 1293 

Wall Temp. Tw 966 943 
Tube Temc) 



Closer to flame 1021 907 888 0.71 133 19 

Closer to wall 979 882 868 0.8 1 111 14 

20 Without other supporting observations (temperature and flow rate of the 

process fluid in the tubes, thermocouples at the target, etc.), it was not possible to 
assess the accuracy of the method ODewitt, 1986). First, the calculated view factor 
from wall to target (F^.J was lower if the target was in closer proximity to the 
burner. This is reasonable because the view factor from flame to target is (l-F^J 
25 by the summation rule. Second, the reflection errors, (Txi-T,), are less at 3.9 

than at 0.9 Mm, which shows that the longer wavelength pyrometer is less sensitive 
to the reflected radiation. The errors also increase if the target is closer to the 
flame in agreement with the concept that errors are mainly caused by the reflected 
component originating from the flame. 

The most questionable deficiency in the above system is the use of T^; as the 
true wall temperature even though it was determined by a pyrometer. The errors 
caused by the reflected radiation from the flame in T^; are comparable to the errors 
in Ts, and therefore should not be ignored. Moreover, the simplifying assumptions 
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(panicularly through the use of a view factor) severely limit this method (Dewitt, 
1986). 

The work done by Suarez-Gonzalez et al. (1987, 1992) also relates, in 
principle, to previous pyrometry approaches. It is the most recent technique 
5 published which discusses measurement of a surface temperarure in the presence of 
hotter surroundings. The problem is again the measurement of the temperature of 
turbine blades in a jet engine. The spectral bands in the two-color system were 
similar to what Atkinson and Strange (1976) used, 0.4 to 1.05 fim, unfiltered, and 
0.4 to 0.85 fitn, filtered. The spectral bands used in their three-color system were a 
10 short wavelength band at 0.3 to 0.85 /im, a middle wavelength band at 0.85 to 1.0 
fim, and a long wavelength band at 1.0 to 1.8 ^m. 

For the two-color system, the unfiltered and filtered signals and Ef were 
produced from the pyrometer targeting the turbine blade. The pyrometer signals, 
and Ef, were the sums of emitted and reflected components. 
15 K-^E^^E^, 



20 



E.^E^^E^ 

The ratio between the filtered and unfdtered signals for a given reflected intensity 
was assumed constant 

I Efr = constant = 

Consequently, the equation for the total unfiltered signal was written as 



(C) 



(D) 



= ^ub + KfEfr (E) 

25 From equation (C), the filtered, reflected signal becomes 

^fr-^f'^Jb (F) 

Combining the preceding two equations results in 

where, if R^f is known, the left-hand side is a function of detectable optical signals 
(including components from both the blade and the flame), while the right-hand side 
is only a function of blade temperature. 
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Since the flame (blackbody) temperature can be estimated on the basis of 
material and energy balances, and R^f can be calculated from equation (D), the term 
(Ey-RyfEf) in the left-hand side of equation (G) is then obtained from the pyrometer 
readings. As a result, the blade temperature can be calculated since a monotonic 
5 relationship exists between (E^b-R^jfEfb), which is equal to (E^-R^^) by equation (G), 
and Tb, within a realistic blade temperature range (< 1650*'C). Here the emissivity 
of the blade and the effect of the engine geometry are considered known. The 
central idea of this approach is to reduce the two equations (C) to one equation (G) 
in order to establish the relationship between the pyrometer readings (their 

10 combination, actually) and the blade temperature, Tb. As a result, only one 
unknown, Tb, can be solved for. 

On the basis of the two-color system, a three-color system also was 
developed. Two combinations, short-middle wavelength bands and middle-long 
wavelength bands were separately used in the same way as equations (C) through 

15 (G). Hence, this three-color system essentially included two, two-color systems and 
therefore allowed two variables, T^, and the flame temperature, Tf, to be measured. 
Unlike being specified based on the estimated flame temperature, as was done for 
in the two-color system, and were now calculated based on a given Tf. 
If Tf was not properly chosen, the resulting incorrect values of R^^, and R^ would 

20 give different blade temperatures by the two systems, which is not realistic. The 
value of Tf (and then R^^ and R^, accordingly) can be chosen iteratively until the 
difference in the Tb's from the two systems becomes acceptably small. 

The most severe limitation of the systems proposed by Suarez-Gonzalez et 
al. is the need for a specified emissivity of the surface. This emissivity is used to 

25 obtain the relationship between the temperature of the blade, T^, and the pyrometer 
signals originating from the blades, Eb. It is also used to obtain the relationship 
between the temperature of the flame, Tf, and the reflected signals, E^^ and Ef,. The 
requirement of a emissivity prevents this technique from being used in many 
applications where the emissivity of the surface is not known. 

30 Measurement of surface temperature has many industrial applications. For 

example, temperature measurement of particles in a furnace may prove very 
beneficial in controlling the operation of the furnace. As an example, the 
temperature of the slumping beds of solids in rotary kilns is a key parameter in 
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hazardous waste incineration. However, direct bed temperature measurement in 

full-scale rotary kilns is difficult due to the rotation of the kiln and the harshness of 

the kiln environment. 

Because of this, an indirea approach is frequently applied in which the 
5 temperature of the kiln combustion gases is measured and used to control the kiln's 

operating conditions. It is assumed that the bed solids' temperature is maintained 

within normal limits, if the proper combustion gas temperature range is maintained. 

This approach is not as attractive as a direct measurement of the solids' 

temperature, especially Lf the solids' temperature can be directly measured a 
10 significant distance before the solids exit the kiln. Direct, upstream measurements 

increase the possibility of actually altering kiln conditions to correct inappropriately 

cool solids temperatures before the solids leave the kiln. 

One possible technique for direct measurement is infrared optical pyrometry. 

At first glance this appears to be a feasible technology if air purging can keep lenses 
15 clean and water cooling can prevent the electronics from overheating (Andreyuk, 

1989); however, an optical method might not obtain the tme bed temperature for the 

following reasons: 

A) The wavelengths selected for optical measurements must avoid 
the principle absorption and emission bands of the major kiln gas 

20 species such as H2O and CO2. 

B) An infrared pyrometer measures the radiation intensity from a 
surface which is the sum of emitted and reflected components. If a 
relatively cool surface is placed in hot surroundings, then the 
reflected component of the intensity will cause the temperature to be 

25 overestimated. 

C) The sight path of the pyrometer might be obscured by small, 
suspended particles such as soot, entrained fines, and smoke, 
therefore the optically determined temperatures will involve errors. 
Conventional pyrometers are not able to correct for the effects of reflected 

30 radiation from the surface of interest. This interference results in apparent surface 
temperatures that are too high. 
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DISCLOSURE OF THE INVENTION 
A method of using multiple single-color pyrometers to correct for the effects 
of reflected radiation when optically measuring the temperature of a relatively cool 
surface in hotter surroundings is disclosed. This technique requires that the 
5 emissivity of the surface be known. 

The multicolor pyrometer measures the intensity of radiation leaving a 
surface over two or more narrow wavelength bands. The intensities permit the 
calculation of the surface's temperature, the effective temperature of the 
surroundings, and in some cases the emissivity of the surface. 
10 The invention differs from the pyrometric techniques previously available in 

its ability to discriminate between reflected and emitted radiation. The method 
obtains measurements of detected radiation at different wavelengths and then by 
solving simultaneous equations relating the unknowns determines the temperature of 
a surface. In the version of the method which utilizes two-color pyrometric 
15 measurements, the emissivity of the surface being measured must be independently 
determined or otherwise known. In the three-color approach, the surface emissivity 
is a third unknown which is solved for along with the surface temperature and the 
effective surrounding's temperature. The three-color method does not require the 
calculation of surface emissivity when measurements at three wavelengths are 
20 utilized asnd the surface is grey over the wavelength range of the measurement. 

Applications include the measurement of the temperature of the clinker 
exiting cement kilns and the temperature of the solids exiting solid hazardous waste 
incinerators. 

In conjunction with the multicolor system, a new method is disclosed to 
25 measure the spectral, surface emissivity of porous materials. This method involves 
heating a sample of the material to a uniform temperature and then suddenly 
exposing it to cool surroundings. Knowing the initial temperature and the optically 
measured brightness temperature enables the calculation of the emissivity. This 
techniques requires that the pyrometer readings be extrapolated back to time zero. 
30 The application of the multicolor techniques to industrial applications, such 

as measuring the temperature of particles, e.g., solids in furnaces and, in particular, 
to incinerators, is complicated by the presence of soot and dust, and by wide 
variations in the types of surfaces present. Experimental results using the two-color 
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approach indicate that interferences from smoke, soot, and dust can have significant 
effects on pyrometer readings. Pyrometers operating at longer wavelengths are less 
sensitive to these interferences. 

While a two-color method is disclosed in extensive detail the methodology of 
5 that method has been extended to a novel, three-color technique which is applicable 
even when the emissivity of the surface is unknown, provided that the emissivity is 
independent of wavelength over the range of wavelengths being used. That is, this 
new method requires that the three spectral surface emissivities are the same. 
Computer simulations show that the three-color system is a promising method. 

10 The three-color technique requires properly selected wavelengths. The 

system is highly sensitive to errors in the individual pyrometer readings if the errors 
are not in the same direction and are not of roughly the same magnitude. The 
three-color system is less sensitive to errors which are in the same direction and of 
approximately the same magnitude. A detailed sensitivity analysis indicates that it 

15 is theoretically favorable to design and fabric a three-color system which uses only 
one lens/ sensor-amplifier-display system for the three-spectral signals. Further . 
analysis suggests that the sensitivity to errors can be reduced 5-20 fold if the signals 
from the instrument are not processed into brightness temperatures, but are used 
directly in the calculation of the temperatures and emissivity. 

20 

BRIEF DESCRIPTION OF THE FIGURES 
FIG. 1 schematically depicts emitted and reflected components of radiation 
detected by an optical pyrometer. 

FIG. 2 schematically depicts a cross section of the 1.8 kW furnace for 
25 emissivity determination showing the oxidi2:ed Cu plate with two thermocouples 
inserted for temperature measurement. 

FIG. 3 schematically depicts a device used to measure the spectral emissivity 
of sand or clay. The insulating cover is quickly removed and the transient 
temperature of the sample is measured optically. 
30 FIG. 4 schematically depicts the top view of the water-cooled sample tray 

with two thermocouples to estimate the copper plate surface temperature. 

FIG. 4 A schematically depicts the same sample tray as FIG. 4 further 
illustrating that the same tray can be used for sand or clay. 

RECTIFIED SHEET (RULE 911) 
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FIG. 4B schematically depicts a side view of the water-cooled sample tray of 
FIG. 4 The figure in the circle shows that the same tray can be used for sand or 
clay. 
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FIG. 5 schematically depicts a rotary kiln simulator. 

FIG. 5A is a cross sectional view of FIG. 5 taken along section lines A-A. 

FIG. 6 is a schematic of optical measurement of spectral, brightness surface 
temperature by two single-color pyrometers. 
5 FIG. 7 is a graph depicting the normal, spectral emissivities of oxidized 

copper plate at 1.05 and 2.2 ^m. 

FIG. 8 is a graph depicting oxidized copper data for two-color pyrometer 
verification. 

FIG. 9 is a graph depicting the normal, spectral emissivity of sand at 1.05 
10 and 2.2 ^m. 

FIG. 10 is a graph depicting the normal, spectral emissivity of Smokey Joe 
clay at 1.05 and 2.0 /xm. 

FIG. 11 is a graph depicting temperature data from water cooled tray test in 
rotary kiln. The thermocouple solids' temperature is measured 1.5 mm below the 
15 surface of the sand. 

FIG, 12 is a graph depicting the optical measurement results at batch, 
transient conditions in the kiln using sand at 2 ipm rotation rate and 5 % fiD 
fraction. 

FIG. 13 is a graph depicting an illustrative example of the spectral, 
20 brightness temperature measurement at the surface temperature 1234*'F. The 

brightness temperatures and the surface temperature then can be used to calculate 
the spectral emissivity. 

FIG. 14 are bar graphs dq>icting sensitivity analysis for errors in pyrometer 
reading for a surface with an emissivity of 0.9. 
25 FIG. 15 are bar graphs depicting sensitivity analysis for errors in pyrometer 

reading for a surface with an emissivity of 0.3. 

FIG. 16 are bar graphs depicting sensitivity analysis for errors in emissivity 
for a surface with an emissivity of 0.9. 

FIG. 17 are bar graphs depicting sensitivity analysis for errors in emissivity 
30 for a surface with an emissivity of 0.3. 

FIG. 18 is a graph depicting the effects of smoke interference to a 1.05 fim 
pyrometer reading for an oxidized stainless steel plate outsider of the kiln at an 
approximate temperature of 1150°F. 

RECTIFIED SHEET (RULE 91) 
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FIG. 18A is a graph depicting the effects of smoke interference to a 2.2 
pyrometer reading for an oxidized stainless steel plate outsider of the kiln at an 
approximate temperature of IISO^F. 
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FIG. 19 is a graph depicting soot interference effects at an acetylene flow 
rate of 11.0 SCFH. 

FIG. 20 is a graph depicting entrained fines inierference^ effects when the 
clay particles are fed at a mass flow rate of 62.5 g/min. 
5 FIG. 21 is contour chart of brightness temperatures, Tbxi and Tbx^i 

corresponding to the surface temperature and the effective surrounding temperature, 
Tsurf and T^^„. The conditions are \ = 1.05 fim, X2 = 2.2 /xm, and e^i =0.3, 
£x2 = OA, The plain numbers and lines are for Tbxi and the bold ones are for Tbx2- 

FIG. 22 is a contour chart of T^^^ as a function of the brightness 
10 temperatures, Tg^i and Tbx2- 

FIG. 22 A is a contour chart of T^^ as a function of the brightness 
temperatures, Tg^i and Tbx2- 

FIG. 23 is a graph which illustrates the effect of errors in the brightness 
temperature at 0.8 ^m on the temperatures of the emissivity. No errors are 
15 associated with the other wavelengths. 

FIG. 23A is a graph which illustrates the effect of errors in the brightness 
temperature at 0.8 /im on the temperatures of the surface and the surroundings. No 
errors are associated with the other wavelengths. 

FIG. 24 is a graph which depicts the effect of errors in the brightness 
20 temperature at 1.6 Atm on the emissivity. No errors are associated with the other 
wavelengths. 

FIG. 24A is a graph which dq)icts the effect of errors in the brightness 
temperature at 1.6 /im on the temperatures of the surface and the surroundings. No 
errors are associated with the other wavelengths. 
25 FIG. 25 is a graph which illustrates the effect of errors in the brightness 

temperature at 2.2 /im on the emissivity. No errors are associated with the other 
wavelengths. 

FIG. 25 A is a graph which illustrates the effect of errors in the brightness 
temperature at 2.2 ^m on the temperatures of the surface and the surroundings. No 
30 errors are associated with the other wavelengths. 

FIG. 26 is a graph which depicts the partial derivative, 9Tsyrf*/5TBXi*, i = l, 
2, 3, as a function of the normalized brightness temperature with errors in the 
brightness temperature of j=0, 1,..., 6. 

RECTIFIED SHEET (RULE 91) 
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FIG. 27 is a graph which depicts the partial derivative, dT^^^^/dT^^^i*, i = l, 
2, 3, as a function of the normalized brightness temperature with errors in the 
brightness temperature of j=0, 1,..., 6. ^ 

FIG. 28 is a graph which depicts the partial derivative, de'^/dT^^*, i = l, 2, 
5 3, as a function of the normalized brightness temperature with errors in the 
brightness temperature of j =0, 1 , . . . , 6. 

FIG. 29 is a graph which illustrates that Type A error in the brightness 
temperature errors in three-color results are relatively insensitive to errors of Type 
A which appear on all three pyrometers with the same magnitude and sign. The 
10 scale of normalized error in the brightness temperature is also shown on the top of 
the group. 

FIG. 29 A is a gr^h illustrating normalized Type A error in brightness 
temperature. 
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FIG. 30 is a graph illustrating the relationship between the derivatives, 
dTBXi*/dRxi* and the normalized radiosity of the three-color system. The small 
values of the derivatives indicate that skipping the conversion of to T^xi decrease 
the sensitivity to errors. 

5 

BEST MODE OF THE INVENTION 
It is now known that heat is a form of energy which is transferable from one 
material to another by either conduction, convection, or radiation. AU materials, at 
aU temperatures, down to near absolute zero, radiate electromagnetic energy. The 
10 wavelengths at which this radiation occurs depends on the nature and temperature of 
the material. The temperature of the material determines both the quantity and the 
wavelengths of energy radiated. As the temperature increases, the wavelengths shift 
to shorter values. In combustion applications, the bulk of the thermal radiation 
occurs in the infrared (2-20 /xm) and near infrared (0.7-2 fim) ranges. Therefore, 
15 the optical devices that measure the intensities of thermal radiation are often 
referred to as infrared radiation thermometers (pyrometers). 

The theory of optical pyrometers is based on the Planck distribution, 

20 This expression relates the amount of energy emitted by a perfect radiator 

(blackbody) per unit wavelength, per unit of time, per unit area, at wavelength \, to 
the temperature of the radiator, T. The values of the radiation constants, Cj and 
are 3.742E10' W*;imVm^ and 1.439E10' fim^K. 

The Planck law reduces to equation (2), called the Wien equation, which is 
25 accurate to within 1 percent when XT is less than 3(XX) finfK: 

E.^^C^k-^e-^''-'' (2) 
This range is applicable for the instant method. Because of internal reflection at the 
surface of a material, the intensity of the radiation passing through a surface, and 
30 thus emitted, is less than that given by equation (2). For a nonblackbody, equation 
(2) becomes 
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where the emissivity which is between O and 1 . equation (3) can also be written 
using a fictitious temperature, T^x called the brightness temperature. This is the 
temperature at which a black body will emit the same amount of energy as a surface 
with emissivity and true temperature T: 



If the emissivity is independent of wavelength (gray) or does not change 
significantly over the selected wavelength range, then equation (3) gives 



which means that the temperature can be determined in spite of the fact that the 
emissivity may not be known. Egx can be calculated using equation (1) or (2), A 
single-color pyrometer can be used to determine the surface brightness temperature, 
Tbx. and then E,, can be calculated using equation (4). The ratio of E>, to E^x, as 
shown in equations (3) and (4), is £x- Th^ combination of one radiation ratio 
thermometer and one single-color pyrometer can, in principle, be used to determine 
the emissivity of the surface. 

Radiation originating from extraneous sources and reflecting from the surface 
of a target increases its spectral radiance and causes errors in radiometric 
temperature measurements. When the magnitude of the reflected radiation is 
comparable to or larger than that of the emitted radiance, errors in optical 
temperature measurements tend to be large. At target temperatures below about 
100 to 150**C, reflected radiation is usually the dominant source of error in 
radiation thermometry (Nutter, 1985). Under certain circumstances, such as in a 
billet-reheating furnace, neglecting the effect of reflected radiation may simply make 
the optical measurements useless. (Hottel, 1986). 

When reflected radiation is significant, the analysis becomes more 
complicated because the radiation detected by a pyrometer will include reflected and 
emitted components. Taking radiation emitted and reflected from the surface into 
consideration, pyrometer readings at wavelengths and are given by 



(4) 




R,^ (r,„^) + (l-£,^)£,,,(r_) 



(6) 
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where T^^^ and T„„. are the actual temperature of the surface and the effective 
temperature of the surroundings. Note that T^^^ is a fictitious temperature which 
accounts for the temperatures of the sources of the reflected radiation and their 
geometric disposition. 

Two wavelengths are necessary to give two equations in the two unknown 
temperatures, T^^^ and T^^. The emissivities in the system (6) are still 
undetermined by these relations and must be obtained from independent 
measurements. The readings R^i have the same units as E^j, and (l-^xi) is, the 
spectral reflectivity of the opaque surface. By Kirchhoffs law, the reflectivity and 
emissivity sum to one. 

Obtaining the emissivity of solid materials is difficult since depends on 
angle, wavelength, and temperature of the surface. One approach to eliminating the 
need to know ex, and ex2 involves taking readings at three wavelengths 
instead of two. The resulting three equations contain five unknowns: 

where, again, T^rf and T^„ are the actual temperature of the surface and the 
effective temperature of the surroundings. 

If the surface can be assumed gray at the wavelengths of interest, i.e., e^x = 
^x2 = fix3 = equations (7) can be solved since there are only three unknowns, 
namely T^, T„„, and e. The system is highly nonlinear and transcendental, and 
must be solved numerically. The assumption of identical spectral emissivities, 
although somewhat restrictive, is reasonable if the wavelengths and the applications 
are carefully selected. In many cases, the emissivity of a surface can be treated as 
a constant if the wavelengths at which the measurements are taken are within a 
narrow range. For example, the spectral, normal emissivity of heavily oxidized 
stainless steel at 1200 K is close to 0.8, over a wavelength range from 0.4 fim to 2 
fim (Incropera and DeWitt, 1985). 

One potential application of pyrometric analysis is the sensing of 
temperatures of solids, particularly solids in closed containers, e.g., solids contained 
in heated kilns. For example, the temperature of the bed of solids in rotary kilns is 



) 

rr / 
r) 
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a key parameter in hazardous waste incineration. However, direct measurement of 
the bed temperature in full-scale rotary kilns is difficult due to the motion and 
structural configuration of the kiln as well as the harshness of the kiln environment. 
Because of this, an indirect approach is frequently applied in which the 
5 temperature of the combustion gases at the exit of the kiln is measured and used to 
control operating conditions. It is then assumed that the temperature of the solids is 
maintained within normal limits if the proper combustion gas temperature range is 
maintained. This common practice is often not reliable. Significant discrepancies 
often exist between the temperature of the exit gas and the temperature of the 
10 solids. 

One possible technique for the direct measurement of the temperature of the 
solids is infrared optical pyrometry. At first glance this appears to be a feasible 
technology if air purging can keep lenses clean and water cooling can prevent the 
electronics from overheating (Andreyuk, 1989). An optical method might not 

15 obtain the true bed temperature for several reasons. The wavelengths selected for 
optical measurements must avoid the principle absorption and emission bands of the 
major combustion gas species such as HjO and COj. An inftared pyrometer 
measures the radiation intensity from a surface which is the sum of emitted and 
reflected components. If a relatively cool surface is placed in hot surroundings, 

20 then the reflected component of the intensity will cause the temperature to be 
overestimated. The sight path of the pyrometer might be obscured by small, 
suspended particles such as smoke, soot, and entrained fines, therefore, the optically 
determined temperatures will affected by particle interference (e.g., emission, 
absorption and scattering) and consequently involve errors, 

25 There exists a need for an optical technique which eliminates interferences 

from hotter surroundings. 

The two-color technique of the instant invention gives rise to a system of 
two equations in four unknowns. 
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These two equations can be solved if the spectral emissivities of the surface are 
known. 

FIG. 1 illustrates that the radiation detected by an optical pyrometer will 
include emitted and reflected components. Here, the surfaces are assumed to be 
5 diffuse emitters and reflectors and the effects of participating gases are neglected. 

To usefully apply equations (6), the proper choices for Xj and must be 
made. If the optical measurement is conducted in the exit region of the Idln, where 
the interferences are not as severe as upstream, and, if the wavelength windows of 
the optical pyrometry are selected such that the emission and absorption bands of 

10 CO2 and H^O can be avoided, the major systematic errors involved in optical 
measurement for solid temperatures wiU then most likely result from reflected 
radiation. Factors affecting these choices include avoidance of the CO2 and HjO 
emission and absorption bands, and the selection of wavelengths that can see 
through clouds of soot, dust, and smoke. 

15 The principle emission and absorption bands for CO2 are 2.64-2.84, 

4.13-4.5, and 13-17 /xm. For HjO, the bands are 2.55-2.84, 5.7-7.6, and 12-25 fim 
(Hottel and Sarofun, 1967). 

The mean penetration distance before scattering is inversely proportional to 
the scanering cross section of the dust, soot, and smoke particles. Since the 

20 diameter of soot particles is small relative to the wavelengths typical of thermal 
radiation, the scanering cross section is proportional to X"^ where X is the 
wavelength (Siegel and Howell, 1981). Increasing the wavelength greatly increases 
the penetration distance before scattering. Hence, longer wavelengths are desirable 
for seeing through clouds of soot. 

25 The fmal factor to consider in selecting wavelengths for equations (6) 

involves the sensitivity of T^^ and T,^^ to errors in the readings, Rxi and R^, By 
taking readings at two wavelengths, equations (6) can be solved simultaneously to 
find both T^^^ and T„„. The design parameters are Xj and Xj. A preliminary 
sensitivity analysis to determine reasonable wavelengths is important. If T^ and 

30 Tsurf, X, and Xj, are specified, equations (6) can be used to calculate RXj and RXj. 
They can then be inverted to obtain T„„ using a numerical solver for a 

nonlinear system. Finally, the inversion process can be repeated with perturbed 
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values of RX, and RXj to check on the sensitivity of the temperatures to errors in 
the readings. Errors of ±5 percent were used in the example that follows. 

The temperatures and wavelengths examined are summarized in Table 2. 
The temperature of the surroundings was set at 1100 K and the temperature of the 
5 surface was set at 900 K. The Newton-Raphson method was used to solve for T^^ 
and Tsu„. Two groups of calculations at wavelengths of 4-5 and 1-2 pLxn were 
performed. These infrared wavelengths were chosen to avoid the principle bands of 
CO2 and H2O. 

Table 2 Summary of sensitivity analysis 



A. Loneer wavelengths (X, 


=4 tim, X,=5 nm. £=0.9) 




Test 


Rm 


Rx2 


Bed Temperature, K 


Solid Temperature, K 


Rm+5% 

Ru-5% 

Standard 

Rv,-5% 


7977 
7217 
7597 
7597 
7597 


5534 
5534 
5534 
5811 
5257 


836.1 
No roots found 

900 

812.3 
No roots found 


1468 
No roots found 

1100 

1496 
No roots found 


B. Shorter waveleneths (X, 


= 1 urn, X,=2 um. e=0.9> 




Test 


Rx, 




Bed Temperature, K 


Solid Temperature, K 


Rm+5% 
Rm-5% 
Standard 
Rx2+5% 
Rx,-5% 


122.1 
110.5 
116.3 
116.3 
116.3 


5242 
5242 
5242 
4980 
5504 


878 
902 
900 
889 
911 


1107 
1093 
1100 
1107 
1091 



25 At Xi = 4 /xm and = 5 ^m, the inversion process yielded two sets of roots, but 
the calculated temperatures were extremely sensitive to small errors in the readings 
Rxi and R^. For the other two cases at 4 and 5 fim the perturbed readings were 
such that no roots were found. Going to shorter wavelengths of 1.0 and 2.0 fim 
gave temperatures which were considerably less sensitive to errors in the pyrometer 

30 readings. For this reason, and because of the availability of instrumentation, 
pyrometers reading at 1.05 and 2.2 ^m were chosen. 

The two-color optical pyrometry technique was tested experimentally in three 
stages. 1) Verification tests were performed using a heavily oxidized copper plate. 
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The plate's emissivity was measured, as described below, and then the same plate 
was placed in a wateijacketed tray and inserted in a pilot-scale rotary kiln at 
elevated temperatures. This permitted the new technique to be tested under 
carefully controlled conditions. Copper was selected for verification purposes 
5 because of its high thermal conductivity. 2) Verification tests were also performed 
using sand and a clay-type soil in the pilot-scale kiln. The tests were conducted for 
both steady state and transient heating conditions. In conjunction with these tests, a 
new method was developed to obtain the required emissivities of the surfaces. 3) 
Interference studies were conducted to simulate the harsh conditions inside actual 

10 rotary kilns. These tests qualitatively examined the effects of smoke, soot, and 
entrained fines on optical pyrometry measurements. 

The optical pyrometer is the key analytical instrument used in this example. 
Two single-color pyrometers were supplied by E^TECHNOLOGY Corporation: the 
PULSAR n. Model 7000SI (1.05 fxm) and Model 7000FF (2.2 ^m). These two 

15 pyrometers operate at wavelengths that avoid the principle emission and adsorption 
bands of CO2 and H2O. 

To measure the emissivity of the copper plate, thermocouples were used with 
the two, single-color pyrometers and a 1.8 kW furnace. FIG. 2 shows a cross 
section of the furnace with the copper plate resting above the heating elements. The 

20 1.8 kW furnace was used to heat a sample to 870°C (1600*^1^ at stead state in 

room-temperature surroundings. In order to eliminate the effects of reflection, the 
sample was heated electrically from below. The two thermocouples measured the 
thermal gradient inside the plate, from which the surface temperature of the 
oxidized copper was estimated by linear extr^lation. Knowing the temperamre of 

25 the surface, T, the desired spectral emissivity was obtained by the procedure 
discussed above. 

When silica sand or a clay-soil were placed in the same apparatus, their 
surfaces were highly nonisothermal, because of their low thermal conductivities, and 
the temperatures of their surfaces were difficult to measure. A new method for 
30 measurement of emissivities was developed using the apparatus shown in FIG. 3. 

The device consists of an electrically heated sample of clay or sand which 
was heavily insulated on the sides and top. A thermocouple inserted at the 
insiilation-sample interface was used to measure the initial surface temperature of 
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the sample after steady state conditions were reached. Then the insulating cover 
was suddenly removed after a pyrometer had been previously focused on the 
surface. The surface temperature was measured optically, with some delay, due to 
the response time ( < 1 second) of the instrument. Meanwhile, a pyrometer, 
5 calibrated against a black body radiator, was sighted on the surface to get the 
brightness temperamre T^, from which was obtained through equation 4. 
According to equation 3, the ratio of E^/Ebx gives eX, which is the desired surface 
emissivity. 

When sand and the clay-soD are placed in the same apparatus, their surfaces 

10 are highly nonisothermal, because of their low thermal conductivities, and their 

surface temperatures are difficult to measure. A new method was developed using 
the apparatus shown in FIG. 3. The device consists of an electrically heated sample 
of clay or sand which is heavily insulated on the sides and top. A thermocouple 
inserted at the insulation-sample interface is used to measure the initial surface 

15 temperature of the sample after steady state conditions have been reached. Then the 
upper layer of insulation is suddenly removed after a pyrometer has been sighted on 
the surface. The surface temperature can be read by the pyrometer aftei* its 1 second 
response time. A data acquisition system was used to record the transient, optically 
determined brightness temperature, T^. The initial surface brightness temperature, 

20 used in determining sample emissivity, must then be obtained by extrapolating to 
time zero. More details are given below. A data acquisition system was used to 
record the transient opticaUy determined brightness temperature, Tg, The initial 
brightness temperature of the surface used for determining emissivity of the sample 
must then be obtained by extrapolating to time zero. 

25 A 130-kW, natural gas-fired, pilot-scale rotary kiln (Owens, 1991) and a 

water-cooled sample tray were employed to simulate the situation in which cool 
solid wastes are placed in a hot environment. A schematic of the tray is given in 
Fig. 4. The tray, with a thermocouple embedded in the sample of sand or clay, 
was kept at a relatively low temperature inside the kiln by the water-cooling system. 

30 The cooled tray can be used for either the copper plate or for other materials with 
minor modifications. The rotary kiln and the associated afterburner are shown in 
FIG. 5. The variable swirl burner can produce temperatures of up to 1620**C in the 
0.61-m-diameter and 0.61-m-long combustion chamber of the kiln. The burner and 
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main section were rotated at speeds up to 2.0 ipm; the exhaust section remained 
stationary. The kiln exhaust gas passed between two opposed, 49 kW natural-gas 
fired burners in an afterburner tower (Owens, 1991). FIG. 6 shews how two, 
single-color pyrometers were used to measure the spectral brightness temperatures 
5 of the sample surface inside the kiln. 

For the interference studies, acetylene (C2H2) was mixed with methane (CH4) 
and burned under fuel rich conditions to generate a sooty environment inside the 
kiln. The afterburner was operated under fuel-lean conditions to ensure a 
sufficiently safe overall fuel/air ratio. In addition, the clay soil, with a particle size 
10 of less than 37 fim, was injected into the kiln through the burner to simulate the 
dusty environment typical of industrial rotary kilns. Commercial smoke emitters 
were used outside of the kiln to examine the interfering effects of smoke in the 
pyrometer sight path. 

An OMEGA® WB-AAI interface card, installed in a Macintosh SE 
15 computer, with the associated software, was used to read and record the output 
from thermocouples and optical pyrometers. The ^jplicable program is attached 
hereto as Appendix. 

FIG. 7 gives the measured spectral emissivities of oxidized copper plate as a 
function of temperature, by the method described above. The difference in the two 
20 thermocouple readings in the plate was typicaUy less than 1**C (the top one always 
read lower as expected). The small temperature gradient implied that the plate was 
close to isothermal condition and that the temperature of the surface could be 
estimated accurately. 

To validate the two-color technique using the copper plate and the water 
25 cooled tray, two single-color pyrometers were used to measure the spectral 

brightness temperatures of the copper plate inside the kiln, as shown in FIG. 6. 
The resulting brightness temperatures were used to calculate the EX's from equation 
(3). The Exi, expressed as R^i, and the corresponding emissivities shown in FIG. 7, 
were substituted in equations (6) and numerical solutions for the temperature of the 
30 surface and the effective temperature of the surroundings were obtained. The 

temperature of the plate was also measured by two thermocouples embedded in the 
copper as shown in FIG. 4. A suction pyrometer was employed to monitor the 
temperature of the gas in the center of the kiln cavity. The wall temperature of the 
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kiln was also continuously recorded during the experiment by a thermocouple 
inserted through the shell of the kiln to within 2 mm of the inner surface. 

FIG. 8 shows the measured temperatures for the kiln and i\}e water cooled 
copper plate inside the kiln. The test lasted for 60-70 minutes and the temperatures 
5 fell by roughly 110°C during the experiment which permitted the optical response 
to temperature variations to be examined. The results are encouraging and show 
good agreement between the optically and directly determined surface temperatures. 
In addition, the optically determined temperature of the surroundings and the gas 
temperatures obtained by the suction pyrometer are in good agreement. 

10 FIGS. 9 and 10 show the measured emissivities of the sand and the clay soil. 

The data are consistent and reproducible. These results were obtained with the 1.8 
kW fiimace. The published spectral emissivity data (Touloukian, 1970) for pure 
crystalline SiOj at 1023**K are 0.4 and 0.063 at 1.0 /xm and 2.0 ^m. The 
discrepancies between the measured and literature values can be explained by the 

15 impurities present in the sand. No published spectral emissivity data for clay-soils 
have been found. 

FIG. 1 1 shows additional results for sand placed in the water-cooled tray in 
the pilot-scale kiln. The temperature of the steel tray was measured by a 
thermocouple soldered to its surface. The temperature of the sand was measured by 

20 another thermocouple located 1.5 mm below its surface so that the optically 
determined temperature of the surface was somewhat higher than the value 
measured with the thermocouple. The optically determined temperature of the 
surroundings is between those of the suction pyrometer and the wall. The optically 
determined temperatures of the sand are quite scattered. A later section will show 

25 that this might be due to the sensitivity of the optically determined temperatures to 
errors in measured emissivities and pyrometer readings, especially when the 
emissivities are low as they are for sand. 

Additional tests were performed with sand by suddenly loading 7.5 kg into 
the kiln. The rotation rate of the kiln was 2 ipm to ensure that the bed was in a 

30 continuous rolling motion and was nearly isothermal. A thermocouple was inserted 
in the sand to monitor its temperature. FIG. 12 shows a comparison of the optically 
measured temperatures of the sand with other temperature in the kUn. The 
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temperature of the sand as measured by the immersed thermocouple and by the 
optical technique are considerably different. 
Emissivitv Measurements 

The accurate measurement of the emissivity of the surface is a key factor 
5 determining the success of the two-color technique. There are some difficulties 
involved in the measurement of emissivity. A conventional radiation ratio 
thermometer might help overcome these difficulties. It could be substituted for the 
thermocouples used in the tray-fiimace studies to get the true temperature of the 
sample's surface, provided that the surface is gray. 

10 However, the error in the temperature of the surface is not expected to be 

significant if a material with a high thermal conductivity is used. Hence, an 
oxidized Cu plate was used for the verification of this novel optical technique. 
Unfortunately, the contaminated solids encountered in thermal treatment applications 
generally have low thermal conductivities. 

15 FIG. 3 shows the apparatus used for the instant method of determining 

emissivity. The instant technique involves suddenly exposing uniformly heated sand 
or clay to room temperature surroundings by rwnoving a layer of insulation. The 
temperature of the rapidly cooling surface is recorded with a single-color pyrometer 
and the initial temp>erature is obtained by extrapolation of the cooling profile to time 

20 zero. FIG. 13 shows an example of this extrapolation scheme for the silica sand. 
Since the ratio of the thermal conductivities between sand and insulation materials 
can range from 5 to 10, the thermal gradients within the sand are less important. 
As a result, readings from a centrally located thermocouple at the sand-insulation 
interface can be approximated as the initial surface temperature, T, before the cover 

25 is removed. This surface T at time zero can be used to calculate E^x by equation 
(1). Immediately after the sample is uncovered, the 1.05 and 2.2 /xm pyrometers 
monitor the dropping spectral, surface brightness temperatures as shown in the FIG. 
13. ITie initial, spectral, brightness temperatures were determined by extrapolating 
these two curves to time zero, Once the three temperatures indicated by arrows in 

30 FIG. 13 were specified, the spectral emissivities of the sand could be calculated in 
the same way as for the copper plate described previously. 

In performing this sudden exposure technique, special attention must be 
given to the following items: a) It takes several hours for the sand to reach a steady 
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state temperature; b) The pyrometer should be focused on the sample's surface 
before the insulation is put in place; c) The output of the thermocouples and of the 
pyrometers should be connected to the same data acquisition system to reduce 
systematic errors. 

5 The analysis of both the water-cooled sand test, shown in FIG. 11, and the 

continuously rolling sand test, shown in FIG. 12, used the emissivities that were 
measured using the sudden exposure technique. Note that the emissivity of the 
surface might change in the rolling bed tests due to the decreased density of the 
tumbling sand. This possible effect has been ignored here. 

10 The results in FIG. 12 show that a large temperature difference exists 

between the surface and the bulk regions of the sand, the surface being much hotter. 
As time increases, this temperature difference decreases. This difference can be 
attributed to the rapid heating of the exposed bed surface by radiation from the 
surroundings during the brief surface exposure interval. 

15 The software which was developed to solve equations 6 (see Appendix) was 

used to investigate the sensitivity of the optically determined temperatures to errors 
in the pyrometer readings or brightness temperatures. FIG. 14 shows typical results 
for a surface with an emissivity of 0.9. For errors in the brightness temperature of 
±10 ®F, the errors in calculated Idln temperatures and solid temperature are 

20 tolerable for wavelengths of 1.05 and 2.2 ^m. 

FIG. 15 shows the same results for an emissivity of 0.3. This value roughly 
corresponds to the values measured for sand. For errors in the brightness 
temperature of ±10 ^'F, the errors in the calculated kiln temperature are still 
acceptable. Unfortunately, the errors in the calculated solid temperature become 

25 significant. This observation may explain the discrepancies seen in the sand data 
shown earlier. These errors can be minimized by careful calibration and operation 
of the instrumentation in the laboratory environment. However, the sensitivity 
shown in FIG. 15 is such that it is doubtful that careful calibration will solve the 
problem. This is especially true in industrial applications where it wiQ be more 

30 difficult to maintain optimum instrument performance. 

In addition to errors in pyrometer readings, the software was used to study 
the sensitivity of the new two color technique to errors in surface emissivity. FIGS. 
16 and 17 illustrate the errors in the resolved kiln and solids temperatures due 
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errors of -±0.02 in emissivities of 0.9 and 0.3, respectively. For an emissivity of 
0.9, the calculated kiln temperatures are more sensitive to these errors, while the 
calculated solids' temperatures are more sensitive at lower emissivities. This can be 
explained as follows. From equations 6 and 7, it is clear that the readings, Rxi, 
5 consist of two parts, one is a function of the solids surface temperature, T^urf, and 
the other of the surroundings temperature, T^^^. For the lower emissivity materials, 
reflection from the surroundings dominates the readings so that the intensity is more 
responsive to the surroundings' temperamre and the optical system becomes a better 
monitor of the surroundings' temperature. On the other hand, a higher emissivity 

10 material makes the pyrometers more responsive to the solids' temperature. 
Transient, Batch Experiments in the PUot-Scale Rotary Kiln 

The results in FIG. 12 show a large temperature different between the 
surface and the bulk temperature of the sand bed, the surface being much hotter. 
As the time increases, this temperature difference decreases. This difference can be 

15 attributed to greater radiation flux from the kiln surroundings initially and/or to 

rapid heating of the exposed bed surface by radiation from the surroundings during 
the brief surface exposure interval. 

A penetration model was used to estimate the surface temperature of the bed 
under these conditions. Radiative transfer is the only heat transfer mode considered 

20 because it dominates at the high temperatures examined. The heat flux to the 

surface, qo", which is assumed constant over the brief exposure interval of less than 
2 seconds, is given by 



25 



30 



qo"^ea(Tj-T,') (8) 

where t is the total emissivity of the surface, <r is the Stefan-Boltzmann constant, 
and T^ and T; are temperatures of the kiln wall and of the surface of the sand at 
time zero. The corresponding temperature of a semi-infinite solid, as a function of 
time and distance from the surface, is given by Incropera and Dewitt (1985) 

At x=0, equation (9) becomes 
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(10) 



where a = k/pc. 

A sample calculation for sand at Tw=870''C and Ti=296''C shows that the 
5 temperature of the surface rises to lOC^C higher than bulk bed temperature in 2 
seconds (See Appendix B) and to about TCC higher in 1 second. Hence, at least 
part of the difference between the temperatures of the bulk and of the surface of the 
bed can be attributed to transient heating of the surface. 
Interi'erence from Smoke. Soot, and Dust 

10 Most of the high temperature experimental work presented above focuses on 

the results from a natural gas-fired, pilot-scale rotary kiln. The gas phase was 
relatively free of particles of smoke, soot, and dust. In industrial-scale rotary kiins, 
pyrometry readings may be compromised by suspended particles. Preliminary 
experiments were performed to look at the sensitivity of the two-color technique to 

15 these interferences. 

The smoke interference tests were conducted outside of the rotary kiln 
simulator using the 1.8 kW furnace with a heavily oxidized stainless steel plate and 
smoke emitters obtained from REGIN HVAC PRODUCTS, Inc. FIG. 18 gives 
results for the two single-color pyrometers at 1.05 and 2.2 /zm. The tests started 

20 with no smoke present, then three smoke emitters placed around the target were lit 
in sequence with approximately three minutes between ignitions. The two 
pyrometers were arranged so that both shared nearly the same sight path. The 
degree of interference varied widely as the smoke plume shifted in and out of the 
sight paths. However, from FIG. 18 it appears that the 2.2 ;im pyrometer was 

25 much less sensitive to the smoke than was the 1 .05 ^m unit. This observation 
supports the conclusions based on scattering theory. 

Incomplete combustion and poor mixing in industrial rotary kilns can cause a 
tremendous amount of soot to form. Optical pyrometer readings are sensitive to 
interferences from soot. An interference study was conducted by injecting acetylene 

30 and natural gas into the kiln at a stoichiometry of 0.7. Fuel lean conditions were 
maintained in the afterburner at all times. The tests were started after the 
water-cooled sample tray with the copper plate was placed in the kiln and the entire 
system was allowed to come to thermal equilibrium. In addition to optical 
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pyrometer readings at 1.05 and 2.2 ^tm, thermocouples inserted in the copper 
yielded temperatures during the experiments. 

A set of tests at acetylene flow rates of 9.2 (5.4), 14.1 (8.3), and 18.7 
(11.0) mVmin (CFH) were performed. The radiation intensities measured by the 
5 pyrometers increased with increasing flow rates. FIG. 19 shows the effects on the 
brightness temperatures and on the optically determined surface temperatures at an 
acetylene flow rate of 18.7 m^/min (11.0 CFH). Temperatures were also recorded 
without any acetylene supplied to the kiln. Introducing the acetylene caused a 
significant increase in the reading of the 1.05 fim pyrometer reading but had little 

10 effect at 2.2 /xm. The higher sensitivity of the 1.05 /xm pyrometer to soot agrees 
with the previous findings for smoke. The optically determined temperatures of the 
surface, based on the measured spectral emissivities for oxidized copper obtained 
earlier (0.9 and 0.85 for 1.05 and 2.2 ^m), decreased as the acetylene flow rates 
and soot levels increased, the maximum decrease was about 33 °C (60°F). To 

15 explore the sensitivity of the measurements to surface emissivity, the optically 

determined temperature of the surface was recalculated with an emissivity of 0.8 at 
both wavelengths. As shown in FIG. 19, the effect of the change in emissivity is 
negligible. The optically determined surface temperature for a low emissivity 
material is expected to be much more sensitive to errors in the pyrometer readings 

20 and to changes in emissivity. Note that the presence of soot increases the brightness 
temperatures at both wavelengths but leads to a reduction in the optically 
determined temperature of the surface. This is counter-intuitive but not surprising 
since the effect of soot (it produces a luminous flame in the kiln) is inteipreted by 
the two-color system as an increase in the temperature of the surroundings. 

25 During the soot study, small amounts of soot deposited on the cool surfaces 

of the sample tray. It is not known if this has a significant impact on the emissivity 
of the copper oxide. 
Interference from Entrained Fines 

The entrainment of fine particles from the burden inside rotary kilns is a 

30 serious problem in hazardous waste incineration and in many other processes. The 
fine particles are transported by a number of different mechanisms including 
susp)ension, saltation, and wall pick-up, and can be carried out of the Idln in the gas 
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stream (Gag, 1991). These suspended particles can interfere with the sight paths of 
pyrometers and can affect their brightness temperature readings. 

In pDot-scale rotary kiln tests, fine clay-soil (< 37 /xm in diameter) was 
added to the air stream of the burner at three different rates, 14.6, 42.1, and 62.5 

5 g/min. This increased the luminosity of the flame. The experimental set-up was 
similar to that for the soot experiments. Before the tests started, the sample tray 
was allowed to come to thermal equilibrium and the dust injection air was left on 
before and after dust addition so as not to change the temperatures in the kiln. The 
effects of the dust on pyrometer brightness temperatures increased with increasing 

10 feed rates. FIG. 20 shows temperatures at a feed rate of 62.5 g/min. The effects 
on the pyrometer readings are significant and the interference at 1.05 yzm is more 
severe than that at 2.2 /xm. The effects of these changes on the optically 
determined surface temperature are small. Using the previously measured 
emissivities (0.9 and 0.85 for 1.05 and 2.2 /xm), the clay caused the optically 

15 determined surface temperatures to increase by only 9*^0 (16*^10, however a 

significant discrepancy exists between the optical and thermocouple readings for the 
temperatures of the copper plate. The interpretation of these results is complicated 
because a thin layer of clay covered the plate's surface and this certainly changed its 
spectral emissivities. For this reason, a set of lower emissivities (0.55 and 0.50 for 

20 1.05 and 2.2 ^m) were used to simulate the effect of the dust layer. As shown in 
FIG. 20, using the lower emissivities gives better agreement with the true 
temperature of the surface but introduces a larger error due to the reflected 
radiation. Note that the variations in T^^ due to the changes in the brightness 
temperatures caused by the dust are in opposite directions if the lower emissivities 

25 are used. 

Three-Color Technique 

As indicated above, equations (6) are unsolvable unless the emissivities of 
the surface are known. In some industrial applications, the materials whose surface 
temperatures are of concern have consistent and unchanging properties over fairly 
30 broad wavelength bands. Their surface emissivities can be determined by prior, 
independent measurements. As a consequence, the two-color method described 
above may be employed to obtain both the true surface temperature and the 
effective temperature of the surroundings. Examples of such surfaces include steel 
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and cement clinker. However, other processes such as hazardous solid waste 
incineration, where the correction of pyrometer readings for reflected radiation is 
also desired, must deal with unspecified materials whose physical properties (e. g., 
emissivities) are unknown and highly variable. 
5 A pyrometer using three wavelengths may find its greatest usefulness where 

the emissivities, although unknown, are independent of wavelength over the range 
of interest. 

If a pyrometer is sited on a small target area in a large enclosure, and if the 
target surface is a diffuse emitter and reflector, then the apparent emissive power of 
10 the surface is 



where T;, is the apparent temperature of the surface at wavelength X, ex is the 
emissivity of the surface at wavelength X, is the actual temperature of the 
surface, is the effective temperature of the surroundings, and E^x is the emissive 
power of a black body at wavelength X. When T, - T^^,, equation 1 gives = Tx. 
15 The quantity, - T„ defines a temperature measurement error. 

Measurements at three wavelengths yield three equations in five unknowns: 

If we also assume tliat the surface is grey then 

e, -€ 

Aj ^2 *3 

and there are now three equations in three unknowns. Solving the three equations 
simultaneously gives the temperature of the surface, T,, t.he emissivity of the 
20 surface, £, and the effective temperature of the surroundings, Tsur 

Measurements at two wavelengths give two equations in two unknowns 
provided that the emissivity of the gray surface is known from an independent 
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source. A two-color approach (at roughly 1 and 2 ftm) was tested in our laboratory 
and was shown to work satisfactorily over a wide range of combustion conditions. 

Sample calculations using the above three equations and selected values of 
wavelength, emissivity, and temperature are shown in Tables A and B. The 
5 wavelengths were chosen with due consideration of CO^ and H^O absorption bands. 
Table A gives three sets of wavelengths, emissivities, and temperatures. Table B 
gives the corresponding apparent temperatures and measurement errors. 

Table A. Three sets of arbitrarily selected conditions. 



Run 1 




1000 




1300 


e 


0.3 


Run 2 


T,(K) 


800 


(K) 


1100 


£ 


0.7 


Run 3 


Ts(K) 


1200 


T„.(K) 


1300 


e 


0.5 
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Table B. The apparent temperatures and temperature measurement errors at 



the conditions of Table A. 



Run 1 






* 


X (jim) 




Apparent Surface 
Temperature 
T, (K) 


Temperature 
Measurement Error 
Tx - T, (K) 


0.6 
1.05 
2.2 


3.281E+1 
5.5168E+3 
3.6552E+5 


1275.5 
1259.5 
1234.9 


275.5 
259.5 
234.9 


Run 2 


0.7 
1.0 
2.1 


5.1218E+0 
2.3785E+2 
6.6497E+3 


1033.5 
1008.5 
947.9 


233.5 
208.5 
147.9 


Run 3 


0.8 
1.3 
2.3 


7.356E+2 
1.5074E+4 
3.9718E+5 


1261.8 
1256.7 
1253.1 


61.8 
56.7 
53.1 



The system of equations (6) can be extended to give three equations in three 
unknowns. 



20 



(7') 



Note that in equations (7'), the previous system with three equations in five 
25 unknowns has been reduced to system where, assuming Rxi, Rx2 and Rxa are known, 
only three variables, e, T.^^ and T^^ are unknown, equations (?') are therefore 
numerically solvable provided that the emissivity of the surface is independent of 
wavelength. The approach embodied in equations (7') has the advantage that if the 
surface changes, an accurate determination of its temperature is still possible. 

30 
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Computer Simulation Results 

Computer simulations of the three-color idea have been performed to gain 
insights into the feasibility of the approach. Again, the Newton-Raphson method 
was used to solve equation (7')- 

Selecting the wavelengths at which the measurements are performed is 
important. The following general rules were taken into consideration. 

1) The principle emission and absorption bands of CO2 and H^O need to 
be avoided; 

2) The assumption of a constant emissivity is better if a narrow range of 
wavelengths is selected; 

3) The spectral intensities based on Plank law of the visible portions of 
thermal radiation decrease rapidly with decreasing wavelength; 

4) The system is more sensitive to experimental errors when longer 
wavelengths are chosen. 

Rule (2) above requires that the three wavelengths are as close together as possible. 

Table 3 shows the sensitivity of the temperatures determined by solving equations 

(7') to the wavelengths selected. 

Table 3 The effect of wavelength selection on the ability 
to accurately invert equations 7 



Wavelength 


Brightness 
Temperature 


Variable 


Correct Value 


Initial Guess 


Resulting Value 


1.6 fim 
2.2 ^im 
2.4 |im 


1415.5 K 
1409.7 K 
1408.5 K 


Tsurf 
Tsun 
Emissivity 


BOOK 
1500 K 
0.5 


1200 K 
1600 K 
0.3 


1341.2 K 
1553.8 K 
0.72 


Case No. 2: 


Wavelength- 


Brightness 
Temperature 


Variable 


Correct Value 


Initial Guess 


Resulting Value 


0.8 Jim 
1.6 |im 
2.2 Jim 


1434.6 K 
1415.5 K 

1409.7 K 


Tsurf 
Tsurr 
Emissivity 


1300 K 
1500 K 
0.5 


1200 K 
1600 K 
0.3 


1293.8 K 
1496.6 K 
0.48 



25 



30 



In this table, results for two sets of wavelengths are presented: 1.6, 2.2, 2.4 
/xm and 0.8, 1.6, 2,2 /xm. The results of the inversion calculations show that the 
second set of wavelengths gives a more accurate solution for T„rf, and T^^^. All 



BNSDOCID: <WO 9704292A1 J_> 



wo 97/04292 



PCT/US96/120O8 



-31- 

of the errors in this calculation results from rounding the brightness temperatures to 
one decimal point. In making this comparison, the same set of initial estimates for 
£, Tsurf, and T^^„ are used. The greater accuracy of the second set is not suiprising 
since, in set No. 1, wavelengths 2,2 and 2.4 /xm are so close together that the 
resulting brightness temperatures, 1409.7 and 1408.5 K, are not sufficiently distinct 
from each other. Hence, Rule (2), which attempts to ensure that the emissivity will 
be constant, can introduce large inaccuracies when using equations (7') due to the 
highly nonlinear nature of the system. 

Using the three wavelengths, 0.8, 1.6 and 2.2 ftm. Table 4 gives computer 
simulation results based on equations (7') for three different sets of temperatures 
and emissivities. 

Table 4 The computer simulation results using the three-color technique 



Brightness 
Temperature 


Variable 


Correct Value 


Initial Guess 


Resulting Value 


Case A: 


1180.1 K 
1174.9 K 

1173.2 K 


Tsurf 
Tsurr 
Emissivity 


HOOK 
1200 K 
0.3 


1400 K 
1600 K 
0.3 


1089.7 K 
1 197.5 K 
0.26 


Case B: 


1451.5 K 
1379.9 K 
1357.1 K 


Tsurf 
Tsurr 
Emissivity 


1200 K 
1600 K 
0.7 


HOOK 

1700 K 
0.5 


1200.2 K 

1600.3 K 
0.70 


CaseC: 


1396.1 K 

1384.2 K 
1380.8 K 


Tsurf 
Tsurr 
Emissivity 


1300 K 
1450 K 
0.5 


1400 K 
1550 K 
0.3 


1300.2 K 
1450 K 
0.50 



Note that, for each case, the three brightness temperatures are calculated from the 
radiosity (emissive and reflective) the pyrometers detect. Each brightness 
temperature is rounded off to the first decimal point. The above simulation results 
clearly demonstrate the potential of this novel three-color technique in measuring the 
correct temperature of surfaces with unspecified emissivities where reflection from 
hotter surroundings cannot be overlooked. 
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Error Analysis 

The above calculations assume that there are no errors in the pyrometer 
readings or their corresponding brightness temperatures, other than those introduced 
by rounding the brightness temperatures to one decimal place. As indicated above, 
5 the calculated temperatures are subject to errors in the pyrometry readings. The 
sensitivity of the optically determined temperatures and emissivities to errors in the 
brightness temperatures is now considered. 

In the above discussion of the two-color system, comparisons between the 
values of the temperatures measured by thermocouple to the calculated temperatures 
10 based on pyrometer readings implied that the resulting temperatures were sensitive 
to experimental errors. In Table 3, the values of Rxi (i=I, 2, and 3) are assumed 
to be exact with no error other than that due to rounding. The results presented in 
Table 3 show that the solution to equations (7') is extremely sensitive to errors in 
the readings, Rxi. 

15 In practice, some degree of experimental error in the optical pyrometer 

readings is inevitable. Given the sensitivity of the solution to these errors, frequent 
or even on-line calibration is desired. A more fundamaital error analysis is now 
presented in order to gain insight into the sensitivity of the multi-color techniques to 
errors. The examples presented here are for a two-color system, merely for the 

20 purpose of illustrative convenience. Similar examples can be prepared for 

three-color systems with little difficulty except that they are difficult to present 
graphically. 

Calculation were performed with the surface emissivities £xi=0-3 at X, = 1.05 
/xm and e^^OA at X2=2.2 /xm. The true surface temperature ranged from 1200 to 

25 1400 K while the effective temperature of the surroundings ranged from 1400 to 
1600 K. FIG. 21 shows a contour chart of brightness temperatures at Xj and Xj. 
With any combination of T^rf and T^„, a corresponding set of two brightness 
temperatures, can be found with this chart. For instance, at Ta,rf= 1270 K and 
T^^=1440 K, the brightness temperatures are, Tbm « 1400 K and Tb;^« 1375 K. 

30 With the same chart in hand, one might conclude that this process of finding Tgxi 
and Tbx2 can be inverted and, given a set of brightness temperatures, the true values 
of T^^rf and T,„„ could be obtained just by looking at FIG. 21. 
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Unfortunately, for most applications, things are more complicated than this. 
FIG. 22 gives a picture of typical inversion calculations and the problems associated 
with them. In order to illustrate those results effectively, two sq)arate graphs are 
given with the temperatures, and T^, represented by shades of gray. The 
5 graphs show that all of the roots of the system, corresponding to the meaningful 
values of T^^rf and T^^„, are found within a narrow band. If there are significant 
errors associated with Tg^i and T^^, the intersection of these two brightness 
temperatures may faU outside this band, and no roots will be found. FIG. 22 also 
provides insight into how the roots respond to errors or changes in the brightness 

10 temperatures. For this particular case, simultaneously positive or negative errors of 
approximately the same magnitude, hereafter referred as Type A errors, may lead to 
small changes in T^^, Likewise for T^, if simultaneously positive or negative 
errors in the brightness temperatures occur, with an approximate ratio of 2:1 for 
ATbx2 to ATbxi, the change in T^„ is small. In contrast, considerable errors can 

15 occur in both temperatures, or no roots may be found at all, as a result of errors in 
Tbxi and Tbx2 ^at are in different directions. Hereafter, errors of this type are 
referred to Type B errors. 

The stability of a three-color system can be treated similarly. The above 
discussion can be repeated with three, three-dimensional graphs, since three 

20 brightness temperatures and one desired variable, either T^^, T^^^ or £, are 

involved. This three-color system will be explored below and it wiU be shown that 
the diameter of the three-dimensional, roughly cylindrical region in which 
convergence can be obtained is small, i.e., the system is sensitive to errors. The 
orientation of the region is diagonal which indicates relative stability if Type A 

25 errors in the T^xi are involved. An analysis of the sensitivity of a multi-color 
system to enors is necessary in order to specify the precision required of the 
proposed instrument. 

Bearing the previously explored concepts in mind, a number of FORTRAN 
programs were developed to solve equations (1) and/or (5). These programs were 

30 used to investigate the sensitivity of the optically determined temperatures and/or 
the surface emissivity to errors in pyrometer readings. 
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Two-Color Technique 

FIG. 14 shows typical results for a surface with an emissivity of 0.9, 
corresponding to the oxidized copper used in the tests described earlier. The 
emissivity is assumed independent of wavelength. For errors in a single brightness 
5 temperature of ±5.6°C (± lO'^F), the errors in the optically determined 
temperatures are within 2 percent for wavelengths of 1.05 and 2,2 fim, 

FIG. 15 shows similar results for an emissivity of 0.3. The emissivity is 
again assumed independent of wavelength. This value roughly corresponds to the 
values measured for sand. For errors in the brightness temperature of ±5.6*C 

10 (±10**F), the errors in the optically determined temperature of the surroundings are 
still acceptable. Unfortunately, the errors in the temperature of the solid become 
significant. This observation may explain the discrepancies seen in the sand data 
shown earlier. These errors can be minimized by more frequent calibration of the 
instrumentation in the laboratory and by on-line calibration in the field. 

15 The sensitivity analysis presented in Fig. 14 and 15 is for Type B errors. 

Errors of this type are expected here since two separate, single-color pyrometers 
were used. 

In addition to errors in pyrometer readings, the software was used to study 
the sensitivity of the two-color technique to errors in the emissivity of the surface. 

20 FIGS. 16 and 17 illustrate the temperatures of the surroundings and of the surface 
due to errors of ±0.02 in baseline emissivities of 0.9 and 0.3, at both wavelengths. 
For an emissivity of 0.9, the optically determined temperatures of the surroundings 
are more sensitive to these errors, while the temperatures of the solids are more 
sensitive at lower emissivities. This can be explained as follows. From equations 

25 (6), it is clear that the readings, R^i, consist of two parts, one is a function of the 
temperature of the surface, T^^, and the other of the temperature of the 
surroundings, T^„. For the lower emissivity materials, reflection from the 
surroundings dominates the readings so that the intensity is more responsive to the 
temperature of the surroundings and the optical system becomes a better monitor of 

30 the surroundings. On the other hand, a higher emissivity material makes the 

pyrometers more responsive to the temperature of the solids. Test results discussed 
above support this explanation. The two-color system is highly sensitive to the 
errors in the emissivity. 
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Three-Color Technique 

To examine the sensitivity of the three-color approach to unavoidable 
random errors, equations (7') were solved numerically. FIGS. 23 to 25 show the 
calculated changes in T,„rf, T^„ and £, corresponding to the small. Type B errors in 
5 a single brightness temperature or reading at wavelengths 0.8, 1.6, and 2.2 fim. 
The errors in the brightness temperatures range from -6 to +8 K. Errors outside 
this narrow range resulted in no roots to the equations (7') being found. FIGS. 23 
to 25 show that the three-color system is extremely sensitive to all of these errors 
such that it is doubtful that even careful calibration will make the system workable 

10 if a combination of three, separate, single-color pyrometers were used. This is 

especially true in industrial applications where it will be more difficult to maintain 
optimum instrument performance. 

Additional calculations on how errors propagate in the three-color system 
were performed in order to provide guidance on reducing the sensitivity of the 

15 system to errors and to minimize the errors. The following discussion examines 
these issues and it also includes an examination of Type A errors. 

equations (7') expresses the relationship between the independent variables, 
T,„rf, T,„^ and and the dependent variables Rxi- The equations are highly 
nonlinear and must be solved numericaUy. To examine the sensitivity of 

20 and e to errors in the readings or brightness temperatures of the pyrometers, 
equations (7') were inverted and expressed as 

T -T (T T T \ ^''^ 

25 £ = £ ('^Bir'^Bxr'^Bx,) 

where the right hand sides are symbolic expressions for the inverse functions. The 
analytical expressions for these functions do not exist but they can be used to 
examine the sensitivity of the system as follows. 
The total differentials of the system can be written as 



30 








^BX, 




dTsx, 
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^BX, 



dT-sA, 
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Numerically, one can approximate the above total differentials and partial 
derivatives using finite difference methods. Taking T^^ as an example, one can set 
dT.„rf = AT,,rf and, dT^JdT^^,-=AT,,JAT^,j (where Tbx2 and Tbx3 are held 
constant). Forward differencing was used here. Theoretically, the total resulting 
5 error in T,^^^, namely AT^^^, is the sum of three components. Each component is 
the product of a coefficient of error propagation, dT,^JdT^y^, and an error in T^^, 

ATh;,. 

Note that the radiosities (lO^W/m^), the brightness or actual temperatures 
(10'' K) and the surface emissivity (10*^ are of very different magnitudes. For the 

10 convenience of this discussion, normalizing all of the variables was performed as 

follows. In any specific case, once the wavelengths are selected, a set of radiosities 
or their corresponding brightness temperatures were determined by equations (7') 
for any given T,^^, T^^^ and £. For example, for the wavelengths 0.8, 1.6 and 2.2 
/xm, T,„^* = 1300 K, T,„^^ = 1500 K, and £° = 0.4, Rx,"=4692.4 W/im^^fxm) or 

15 Tbx,'*^ 1450.3 K, Rx2- = 67569.0 W/Cm^^^m) or Tbx2^ = 1434.1 K and R;^°=75443.0 
W/(m-»/im) or Tb^'' = 1429.0 K. The symbol C) here stands for the correct value 
of a variable. All of the aforementioned correct values will be used as the 
characteristic values of the corresponding variables in the following discussion and 
the normalized variables will hereafter be denoted by the symbol (*)- For example, 

20 if T.„rt = 1310 K, then T„rf'=T^rf/T,„rf" =(1310 K)/(1300 K) = 1.008; and if 

Tbx, = 1448 K, then Tbxi'=Tm/Tbx,*' = (1448 K)/(1450.3 K)=:0,9984. Hence, the 
normalized variables are all of order 10° and the new, normalized system become 



25 









dTl. 




(13) 


dc' . 







30 The normalized, partial derivatives appearing in equations (13) are plotted in 

FIGS. 26 to 28 as a function of the normalized brightness temperature. For 
example, FIG. 26 gives dT^^//dT^^' as a function of Texi' which is (TBx,'*±j)/TBXi*'> 
i = l, 2, 3 and j=0, 1,...6. The error in the brightness temperatures, given by the 
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integer j, was varied over a limited range because large variations made it 
impossible to solve equations (7'), i.e., it was outside of the range over which the 
roots of the system exist. FIGS. 27 and 28 give similar plots for dT^^„ydT^'' and 

5 Significant insights can be gained from KEGS. 26 to 28. First, the partial 

derivatives are as large as +1000 and -1400. This means that small errors in the 
brightness temperature will be amplified dramatically into the total error estimated 
by equations (13). The severity of the total error depends on the relative location of 
the errors as given by the normalized brightness temperatures. Second, the figures 

10 show that the partial derivatives with respect to T^^' are always negative while the 
other derivatives are always positive. As shown below, this can lead to cancellation 
of errors and improved stability. Third, the error, if it occurs in the single, shortest 
wavelength (0.8 /im), will be relatively unimportant since its coefficients of error 
propagation are smaller than the others. Recall the above discussion, that longer 

15 wavelengths tend to increase the sensitivity of a system to errors and recall Rule 4 
above. Finally, FIGS. 26-28 reveal the reason why errors of Type A are less 
detrimental. Recall that errors of Type A are all in the same direction and are of 
approximately the same magnitude. Errors of Type B are in different directions and 
include errors which affect only one brightness temperature, the other temperatures 

20 being correct. For errors of Type A, the individual errors tend to cancel each other 
with the result that there is a significantly smaller affect on the whole system. On 
the other hand, errors of Type B can lead to large total errors and instabilities. 

The above discussion is critically important to the design of the 
instrumentation. If a three-color system is designed in a fashion similar to the two, 

25 single-color system tested in this study, with three individual units and their separate 
lens/sensor-amplifier display systems, errors of Type B will be unavoidable and the 
system will be extremely difficult, if not impossible, to use. However, if the 
instrument is designed and fabricated in such a way that only one 
lens/sensor-amplifier-display system is shared by all three wavelengths with selected 

30 spectral signals produced from a rotating filter wheel, the stability and sensitivity of 
the system will be much improved. This is because if the frequency at which the 
filter wheel completes a circle is sufficiently high compared to the fluctuation 
frequency of the system, it is reasonable to assume that few errors of Type B will 
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occur if the system is carefully calibrated. The effects of the magnitude of the error 
on the optically determined temperatures and emissivity, as they are most likely 
caused by the fluctuating zero shifts of the lens/ sensor-amp lifier-display system, can 
also be reduced. This can be accomplished if the time-averaged values of the 
5 optically determined temperatures and emissivity are taken over the period 
(reciprocal of frequency) of the system's fluctuation. 

Examples of the calculated insensitivity of the system to errors of Type A 
are shown in FIG, 29 where, as before, the baseline temperatures and emissivity are 
T,„rf=1300 K, T„„=1500 K and £=0.4 at wavelengths of 0.8, 1.6 and 2.2 /xm. 

10 The errors are applied equally to all of the brightness temperatures. FIG, 29 shows 
that the resulting calculated values of T^^.^, T,^„ and s are relatively close to the 
baseline values. The computer simulation results show that the three-color system 
is fairly stable within a range of brightness temperature of ± 80 K. These results 
shows that the three-color technique is a feasible approach if the errors in the 

15 individual brightness temperatures are all in the same direction and of similar 
magnitude. 

In order to further reduce the introduction of errors into the analysis, the 
pyrometer signals at the individual wavelengths, Rxi, need not be converted to 
brightness temperatures, Texi- They can be retained as voltage or amperage signals 
20 which are linearly and directly related to the radiosity. The advantages of doing 
this are shown as follows. 

Plank law is, according to the theory, 

=qA-'(^^^'^^'-i)'' (1) 

25 

equation (I) is for nonblack surfaces since the temperature, Tg, is the brightness 
temperature. If the brightness temperature is also accounting for reflected radiation, 

is replaced with a reading, R^i- The expression of Tsxi as a function of Rxi is 
given by 



wo 97/04292 



PCTAJS96/12008 



-39- 



Normalizing the Tg^j and using the T^^^" and given previously, equation 
(14) can be normalized as 



T = 



+ 1 



Taking the derivative of Texj' with the respect to R^j' gives 



dR 



Xi T*" 0 6 D* 



+ 1 



In' 



+ 1 



10 Since the term dT^yJ" in equations (13) can be written as 



(15) 



(16) 



(17) 



15 



20 



25 



30 



substituting for the dT^iJ" with the above expression, equations (13) become 



9t: 



surf 



di?: 



dr =-^i=^d/?: +^Vd/j: +-^^d/?; 



(18) 



de- =-^d/?: +-^d;?; +^d/?; 



FIG. 30 shows the relationship between the dTex^'/dRxi" and the R^i' 
described by equation (16). The conditions for these calculations are the same as 
those associated with FIGS. 26 to 28. The small value of dTgxiVdRxi' over a wide 
range of R^.* indicates that directly using Rxi as the independent variables, instead of 
the Tbx.; will reduce the sensitivity of the system to errors. The overall error can 
also be minimized since electronically converting Rxj to T^xi and vice versa 
introduces additional errors. The numerical significance of the dTgxiVdRxi* can also 
be determine>d quantitatively since the new coefficients of error propagation, 
dT,^^'/dR^\ aT,„„7aRxi' and de'/dR^', are the products of the old ones, namely 
dT,^//dT^^\ aT^y^VaTexi* and de'/dT^y^*, shown in equations (13), and their 
corresponding dTexiVdRxi*. The data given in FIG. 30 indicate that, at 0.8 and 2.2 
^m, the new coefficients of error propagation will reduce to approximately 
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one-tenth and one-twentieth relative to the previous values. At 1.6 fim, the new 
coefficient will be reduced about five-fold. 

i) Using sand and heavily oxidized copper plate, the feasibility of using 
two, single color pyrometers to correct for the effects of reflected radiation when 

5 optically measuring the temperature of a relatively cool surface in hotter 

suiTOundings has been successfully demonstrated. This new technique requires that 
the emissivity of the surface be known. 

ii) The surface temperature, T^^^, is sensitive to errors in the pyrometer 
readings, especially when the spectral surface emissivities are low. This sensitivity 

10 increases with increasing wavelength. The penetration depth before scattering 
increases with increasing wavelength. Hence, a compromise was reached and 
pyrometers operating at 1.05 and 2.2 fitn were selected. 

iii) A new technique has been developed to measure the spectral, surface 
emissivity for porous materials. This technique involves heating a sample of the 

15 material to a uniform temperature and then suddenly exposing it to cool 

surroundings. Knowing the initial temperature, and optically recording the 
brightness temperature, enables the calculation of the emissivity. This techniques 
requires that the pyrometer reading be extrapolated back to time zero. 

iv) Calculations performed with software developed as part of this project 
20 show that for materials with lower emissivities (<0.4), considerable errors in the 

surface temperatures can occur when applying the new two-color technique. 
Careful calibration and operation can help reduce these errors but their magnitude 
casts some doubt over the practical utility of this technology for low emissivity 
materials in industrial settings. 

25 v) The results of a qualitative, interference study using smoke, soot, and 

dust imply that aU of these materials can have significant effects on pyrometer 
readings. However, the interference study for soot and entrained fines has shown 
that these effects weakly affect the optically determined surface temperature for high 
emissivity materials. For materials with low surface emissivities, the effects of soot 

30 and smoke are expected to be much more severe and possibly even prevent 

converged solutions to equation 6 from being obtained. The experiments indicate 
that the 2.2 fim pyrometer is less sensitive than the 1.05 ftm pyrometer to 
interferences from smoke, soot, and dust. 
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Appendix A: 
The FORTRAN Computer Code 



^23456789 012345 6789 012 345 678901234 56789012345 678901234 5 67890123456789 012 

I AUTHOR: GEOFF SILCOX 
* 3290 MSB 

* DEPARTMENT OF CHEMICAL ENGINEERING 

* UNIVERSITY OF UTAH 

* SALT LAKE CITY, UTAH 84112 

* PHONE: (801) 581-8820 



LAST UPDATED: 1 APRIL 1991 



* THIS PROGRAM ESTIMATES THE TEMPSRATTOES OF SURFACES BY TWO COLOR 
* PYROMETRY TWO OPTIONS ARE AVAILABLE: IFLAG = 1 IMPLIES THAT 

* THE SOIL TEMPERATURE (TS) AND THE EFFECTIVE FURNACE TEMPERATURE 

* (TF) ARE KNOWN AND WE WISH TO CALCULATED READING 1 (Rl) AND 

* [ ^READING 2 <R2) . IFLAG = 2 IMPLIES THAT Rl AND R2 ARE KNOWN AND 
^ . . . . .THAT WE WISH TO CALCULATE TS AND TF . 



REQUIRED INPUT 



IFLAG: IF IFLAG = 1 THEN CALCULATE Rl AND R2 . 

IF IFLAG = 2 THEN CALCULATE TS AND TF. 
WLl = WAVELENGTH 1, MICRONS. 
WL2 = WAVELENGTH 2, MICRONS. 
EMIS = EMISSIVITY OF SOIL. 

TS = TEMPERATURE OF SOIL OR INTIAL GUESS THEREOF, K. 

TF ^ EFFECTIVE FURNACE TEMPERATURE OR INTIAL GUESS THEREOF, K. 

Rl = READING 1 AT WAVELENGTH 1, W/M"2 -MICRON. 

R2 = READING 2 AT WAVELENGTH 2, W/M"2 -MICRON . ^ 
EPS = CONVERGENCE CRITERION FOR NEWTON RAPHSON METHOD, W/M 2. 
IMAX = MAXIMUM ALLOWED NUMBER OF ITERATIONS. 
DELTA = RELAXATION FACTOR. 

THE UNITS ON CI AND C2 ARE W -MICRONS "4 /M" 2 AND MICRONS-K, 
RESPECTIVELY . 



PROGRAM TWOCOL 

OPEN (14, FILES 'TC2. DAT' , STATUS=' OLD' ) 
OPEN(18,FILE='TC2.RES' , STATUS=' NEW' ) 

Cl=3 .742E8 

C2=1.439E4 
READ (14, *) IFLAG,WL1,WL2,EMIS1,EMIS2 

WRITE(18,11) IFLAG, WLl, WL2,EMIS1,EMIS2 

11 FORMAT (XX, 4X/5X, 'IFLAG = ',I2/5X,'WL1 = ',F12,4,' MICRONS 
C /SX.'WL2 = ',F12.4,' MICRONS ' /5X EMIS 1 = ' , E14 . 5/5X, 
C 'EMIS2 = ',E14.5) 

IF (IFLAG. EQ.l) THEN 
READ (14,*) *TS,TF 
WRIT E (18 12) TS , TF 

12 FORMAT (IX, 4X/5X, 'TS = ' , E14 . 5 , ' K'/5X,'TF = ',E14.5,' K' ) 
Rl = RAD (EMIS1,TS,TF, WL1,C1,C2) 

R2 = RAD(EMIS2,TS,TF,WL2,C1,C2) 
WRITE (18 , 13 ) Rl , R-2 

13 FORMAT (IX, 4X/5X, 'Rl = ',E14.5,' W/M"2 - MICRON ' /5X ,' R2 = 
C E14.5,' W/M" 2 -MICRON' ) 



wo 97/04292 



PCT/US96/12008 



•43- 

STOP 

EI^E IF (IFLAG.EQ.2) THEN 

READ(14, *) TS,TF, R1,R2, EPS, IMAX, DELTA 

WRITE (18,14) TS , TF, Rl , R2 , EPS , IMAX, DELTA 
5 14 FORMAT (IX, 4X/5X, 'TS = ',B14.5,' K'/5X,'TF = ',E14.5,' K' 

C /5X, 'Rl = ',E14.S,' W/M" 2 -MICRON' /5X, 'R2 = ',E14.5, 

C ' W/M"2 -MICRON' /5X, 'EPS = ',E14.5,' W/M^2 ' 

C /5X,'IMAX = ' ,I6/5X, 'DELTA = ',F12.4) 

Bl = R1-RAD(EMIS1,TS,TF,WL1,C1,C2} 
10 B2 = R2-RAD(EMIS2,TS,TF,WL2,C1,C2) 

DO 20 I=1,IMAX 
OTS = TS 
CTF = TF 

All =r (1. -EMIS1)*DEDT(TF, WLl, CI, C2) 
15 A12 = EMIS1*DEDT(TS, VrLl,Cl,C2) 

A21 = (l.-EMIS2)*DEDT(TF,WL2,Cl,C2) 

A22 = EMIS2*DEDT(TS,WL2,C1,C2) 

DETJ = A11*A22-A21*AJ.2 

IF (DETJ.KQ.O.) THEN 
20 WRITE (18,23) I 

23 FORMAT (IX, 4X, 'JACOB XAN IS SINGULAR AND I = ',14) 

STOP 

ENDIF 

DF = (B1*A22-B2*A12) /DETJ 
25 DS = (B2*A11-B1*A21) /DETJ 

TF = OTF+DELTA*Dr 
TS = OTS-i-DELTA*DS 

Bl = R1-RAD(EMIS1,TS,TF, WLl, C1,C2) 
B2 = R2-RAD(EMIS2,TS,TF, WL2,C1,C2) 
30 TEST = SQRT(B1*B1+B2*B2} 

IF (TEST.LE.EPS) GO TO 99 
20 CONTINUE 
ELSE 

WRITE (18,21) 

35 21 FORMAT (IX, 4X, ' INAPPROPRIATE VALUE OF IFLAG.') 

STOP 
ENDIF 

WRITE(1S,1S) 

15 FORMAT(lX,4X, ' IMAX EXCEEDED WITHOUT CONVERGENCE') 

40 STOP 

99 CONTINUE 

WRITE (18, 16) TF,TS,I 

16 FORMAT (IX, 4X/5X, 'TF = ',F12.4,' K' /5X, 'TS = ',F12.4,' K' 
C /5X, 'I = ' ,14) 

45 STOP 
END 



FUNCTION RAD(EMIS,TS,TF,WL,C1,C2) 
50 R AD = EMIS*EBL(TS,WL,C1,C2) 

RETURN 
END 



FUNCTION EBL{T, WL,C1,C2) 
55 EBL = CI/ (WL**5* (EXP(C2/ (WL*T) ) -1 . ) ) 

RETURN 
END 
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FONCTION DEDT(T, WL,C1,C2) 

DEDT = C1*C2*EXP (C2/ (WI.*T) ) /(WL**6*T*T* (EXP(C2/ (WL*T) ) -1 . ) **2) 

RETURN 

END 



*234 56789 01234 567 8901234567890123456789 012345678901234 5 67890123456789012 
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Appendix B: 

Sample Calculation of Initial Surface Temperature 

5 Assuming that radiation is the only heat transfer mode imp>ortant here and that 

the bed can be approximated as a semi-infinit solid, the surface temperature of the 
bed can be described by: 

k 

where T(0,t) is the bed surface temperature at x=0 at time t, is the initial bed 
temperature, qo" is the radiative heat flux, and k and a are thermal conductivity and 
10 thermal diffusivity, respectively. 

Assuming qo" is constant over a brief exposure time, it can be expressed as: 

where is the kiln wall temperature, £ is the total emissivity of the bed, and <r is 
Stefan-Boltzmarm constant. 

Taking T^^IASO'^F, Ti=600**F and the published physical properties data for 
15 sand at 300 ^'K (p = 1515 kg/m^ k=0.27 W/m*K, Cp=800 J/kg^K), one can easUy 

calculate that qo" is about 60 kW/m^ AT=118.4 (t^l sec), and AT=167.4 °F 
(t=2 sec). 

The present invention may be embodied in other specific forms without 
departing from its spirit or essential characteristics. The described embodiments are 
20 to be considered in all respects only as illustrative and not restrictive. The scope of 

the invention is, therefore, indicated by the appended claims rather than by the 
foregoing description. All changes which come within the meaning and range of 
equivalency of the claim are to be embraced within their scope. 
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TWO-COLOR SIMULATION FORTRAN CODE 



The following is a listing of computer code for the two-color pyrometry 
simulation. Samples of input and output files are also given. 



( 1 ) Source code of "TWOCOL.FOR" 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C THIS PROGRAM ESTIMATES THE TEMPERATURES OF SURFACES AIvTD SURROUNDINGS 

C BY TWO COLOR PYROMETRY. TWO OPTIONS ARE AVAXUABLE: 

C 

C IFLAG = 1 IMPLIES THAT 

C THE SURPACE TEMPERATURE (TS) AND THE EFFECTIVE SURROUNDING TEMPERATURE 

C (TF) ARE KNOWN. READING 1 (Rl) AND READING 2 (R2J AND THE CORRESPONDING 

C BRIGHTNESS TEMPERATURES (BRTl) AND (BRT2) ARE CALCUALTED. 
C 

C IFLAG = 2 IMPLIES THAT 

C Rl AND R2 (OR, BRTl AND BRT2 > ARE KNOWN. THE TS AND TF ARE 

C CALCULATED. 

C 

c 

C WRITTEN INITIALLY BY GEOFF SILCOX AND REVISED BY DACHUN GAO 

C TO INCLUDE BRIGHTNESS TEMPERATURE CALCULATION AND THE 

C SEQUENTIAL SCREEN TEST FOR SYSTEM SENSITIVITY TO ERRORS. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
*2345678901234567890123456789012345678901234567890123456789012 3456789012 

* REQUIRED INPUT 



.IFLAG: IF IFLAG = 1 THEN CALCULATE Rl , R2 . BRTl, AND BRT2 . 

TF IFLAG = 2 THEN CALCUALTE TS AND TF . 
.WLl = WAVELENGTH 1, MICRONS. 
.WL2 = WAVELENGTH 2, MICRONS. 
.EMISl = EMISSIVITY OF SURFACE AT WLl. 
.EMIS2 = EMISSIVITY OF SURFACE AT WL2 . 
.TS = TEMPERATURE OF SURFACE, K. 
,TF = EFFECTIVE SURROUNDING TEMPERATURE, K. 

Rl = READING 1 AT WAVELENGTH 1, W/M"2 -MICRON . 
.R2 = READING 2 AT WAVELENGTH 2, W/M"2 ^MICRON . 

BRTl = BRIGHTNESS TEMPERATURE FOR Rl , K, 

BRT2 = BRIGHTNESS TEMPERATURE FOR R2 , K. 

ERR = ERROR RESULTING FROM OTHER SOURCES IN TERM OF BRIGHTNESS TEMPERATURE, K. 
EPS = CONVERGENCE CRITERION FOR NEWTON RAPHSON METHOD, W/M*2 . 
IMAX = MAXIMUM ALLOWED NUMBER OF ITERATIONS . 
DELTA = RELAXATION FACTOR. 



THE UNITS ON CI AND C2 ARE W* MICRONS -"4 /M" 2 AND MICRONS- K, RESPECTIVELY. 
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22i 

* 23 4 5 678901234567890123456789012 34 56789012 34 567 890 12 34 56789012345 S789012 
PRCXSRAM TWOCOL 

DIMENSION ERR(lOl) , EMISSl (101) ,EMISS2 (101) ,BRT1 (101) .BRT2 (101) 
OPEN ( 14 , FILE= ' TCI 1 . DAT ' , STATUS= ' OLD * ) 
OPEN{18,FILE=*TCll.RES* .STATUS='NEW' ) 

Cl=3 .742E8 

C2=1.439E4 

READ(14,') N, IFLAG, WL1,WL2 

WRITEC , 11) IFLAG, WL1,WL2 

WRITECIS, 11) N,IFUAG,WL1,WL2 

11 FORMAT (IX, 4X/5X, -N = M2/5X, 'IFLAG = M2/5X.-WL1 = •,F12.4, 
C * MICRONS • /5X, •WL2 = ',F12.4,* MICRONS') 

IF (IFLAG.EQ.l) THEN 

READ (14,*) EMISl , EMIS2 , TS , TF 

WRITE (18,12) EMISl , EMIS2 , TS , TF 
WRITE(*, 12) EMISl, EMIS2,TS,TF 

12 FORMAT ( IX, 4X/5X, 'EMISl = * , F14 . 2 / 5X , ' EMIS2 = •,F14.2 
C /5X,'TS = ',F14.1,' KV5X,-TF = ',F14.1,' K*) 

Rl = RAD(EMIS1,TS,TF,WL,1,C1,C2) 
R2 = RAD(EMIS2,TS,TF,WL2,C1,C2) 
WRITE(18.13) R1,R2 
WRITE ( * 13) Rl R2 

13 FORMAT ( IX. 4X/5X, -Rl '= ',£14. 5,* W/M'^ 2 •MICRON' /5X, ' R2 = 
C E14.5,' W/M*2 •MICRON • ) 

BRTll = BRT(R1, WL1,C1,C2) 
BRT12 = BRT(R2,WL2,C1,C2) 
WRITE (18, 17) BRTll, BRT12 
WRITEC,17) BRTll, BRT12 
17 FORMAT ( IX, 4X/5X, -BRTl = ',F14.1,' K'/5X,'BRT2 = 

C F14.1, • K" ) 

GO TO 155 
ELSE IF (IFLAG.EQ.2) THEN 

READ(14,*) DUMMY, TSCTFO, EPS, IMAX, DELTA 
WRITEdS, 14) TSO.TFO, EPS, IMAX, DELTA 
WRITE (-,14) TSO , TFO , EPS , IMAX , DELTA 

14 FORMATdX, 4X/5X, -TSO = ',F14.1,' KV5X, 

C 'TFO = '.F14.1.' K'/5X,'E?S- ',£14. 5,' W/M''2' 

C /SX, 'IMAX = • ,I6/5X, 'DELTA = •,F12.4) 

DO 111 L=1,N 

READ(14,*) ERR(L) , EMISSl (L) , EMISS2 (L) ,BRT1(L) ,BRT2(L) 
111 CONTINUE 

WRITEdS, 16) 
WRITE {* , 16) 

16 FORMAT {4X, 'ERR* , 6X, 'BRTl ' , 6X, •BRT2' , 6X, 'EMISl ' , 6X, 

C •EMXS2' , 6X, 'TS* ,6X, 'TF* , lOX, ' I* ) 

DO 500 K=1.N 

EMIS1=EMIS$1(K) 
EMIS2=EMISS2 (K) 
TS=TSO 
Tr=TFO 

. Rl = EBL(BRT1 (K) ,WH,C1,C2) 
R2 = EBL(BRT2(K) ,WL2,C1,C2) 
Bl = R1-RA0(EMIS1,TS.TF,WL1,C1,C2) 
B2 = R2-RAD(EMIS2,TS,TF,WL2,C1,C2) 
DO 20 1 = 1, IMAX 

OTS = TS 
OTr = TF 

All = d. -EMISl) •DEDT(TF,Vsl.l, CI, C2) 

A12 = EMISl *DEDT(TS,WL1, CI, C2) 

A21 = d.-EMIS2) •DEDT(TF,WL2,C1,C2) 

A22 = EMIS2»DEOT(TS,WL2,Cl,C2) 

DETJ = A11*A22-A21*A12 

IF (DETJ.EO.O.) THEN 

WRITEdS, 23) K,I 
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WRITEC.23) K,I 
23 FORK.\T (IX, 4X. • JACOBIAN IS SINGUIAR AND ( K, I ) = 

C C • , 14, • , • .14. • ) ' ) 

GO TO 100 
ENDIF 

DF = {B1*A22-B2*A12) /DETJ 
DS = (B2*A11-B1*A21)/DETJ 
TF = OTF+ DELTA *DF 
TS = OTS+DELTA-DS 

Bl = R1-RAD(EMIS1,TS,TF,WL1,C1,C2) 
B2 = R2-RAD(EMIS2,TS,TF,WL2.C1.C2) 
TEST = SQRT(B1*B1+B2*B2) 

IF (TEST.LE.EPS) GO TO 99 

20 CONTINUE 
ELSE 

WRITE (•,21) 
WRITE(18,21) 

21 FORMAT ( IX, 4X, 'INAPPROPRIATE VALUE OF IFLAG . ' ) 
GO TO 100 

ENDIF 

WRITE(18,15) K,I 
WRITE(*,15) K,I 
15 FORMAT {IX, 4X, 'IHAX EXCEEDED WITHOUT CONVERGENCE AND 
C CK,I) = ( M4, ',14, ') •) 
GO TO 100 
99 CONTINUE 

WRITE (*, 19) ERR(K) ,BRT1(K) , BRT2 (K) , EMISl , EMIS2 , TS , TF, I 
WRITE(18.19) ERR{K) ,BRT1(K) ,BRT2{K) , EMISl , EMIS2 , TS , TF , I 
19 FORMAT { IX. F8 . 2 , 2X. F8 . 1 , 2X, F8 . 1 , 2X, F6 . 4, 2X, F6 - 4 , 2X, F8 . 1 , 
C 2X,F8.1,2X,I6) 
100 CONTINUE 
500 CONTINUE 
155 CONTINUE 
CLOSE (18) 
PAUSE 
STOP 
END 

*2345678901234567890123456789012345678901234S678901234567890123456789012 



FUNCTION RAD(EMIS,TS,TF,WL,C1,C2) 

RAD = (1.-EMIS)*EBL(TF,WL,C1.C2)+EMIS*EBL(TS,WL,C1,C2) 

RETURN 

END 

FUNCTION EBL(T,WL.C1,C2) 

EBL = CI/ (WL**5* (EXP<C2/ (WL*T) ) -1. ) ) 

RETURN 

END 

FUNCTION BRT(R,WL,C1.C2) 

BRT = C2/(WL«AL(X;((C1/(R*WL**5)+1.))) 

RETURN 

END 

FUNCTION DEDT(T,VJL.C1,C2) 

DEDT = C1*C2»EXP(C2/(WL-T) )/(WL**6*T*TMEXP(C2/(WL*T) )-l-)**2) 

RETURN 
END 



•23456789012 34 567890123456789012345678901234567890123456789012 3456789012 
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/\ oaliipJC iJUipUL XiiC. 


















N = 11 




















IFLAG = 


2 


















WLl = 


1.0500 MICRONS 


















WL2 = 


2.2000 MICRONS 


















TSO = 


1200.0 K 


















TFO = 


1600 . 0 K 


















EPS = 


0 . lOOOOE+00 W/M*2 


















XMAX = 


1000 




















^ ft ft ft 
. 1000 


















ERR 


BRTl BRT2 


EMI SI 


EMIS2 




TS 




TF 




I 




1 4^ D . 0 13 94 . 7 


5000 


. 5000 




IOTA 

12 30 


- 9 


1517 


. 2 


113 


— 12 . 00 


1426 . 0 13 97 . 7 


5000 


. 5000 




1244 


. 2 


1514 


. 6 


114 


O ft A 


1426.0 1400.7 


5000 


.5000 




1257 


. 7 


1511 


. 6 


114 




14^o.U i.4UJ./ 


5000 


C ftrt ft 

. 5000 




1271 


c 


1508 


. 3 


114 


-3,00 


1426.0 1406.7 


5000 


.5000 




1285. 


.5 


1504 


.4 


114 


.00 


1426.0 1409.7 


5000 


.5000 




1300. 


.0 


1499 


,9 


112 


3.00 


1426.0 1412.7 


5000 


.5000 




1315, 


.1 


1494 


.7 


112 


6.00 


-1426.0 1415.7 


5000 


.5000 




1331. 


.2 


1488 


.3 


113 


9.00 


1426.0 1418-7 


5000 


.5000 




1348. 


.8 


1480 


. 4 


111 


12.00 


1426.0 1421.7 


5000 


.5000 




1369. 


1 


1469 


.5 


110 


15.00 


1426.0 1424.7 


5000 


.5000 




1396. 


5 


1451 


.5 


114 
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APPENDIX D 



THREE-COLOR SIMULATION FORTRAN CODE 



The following is a listing of the computer code for the three-color pyrometry 
simulation. Samples of input and output files are also given. 



(I) Source code of *'TRICOL.FOR" 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C THIS PROGRAM ESTIMATES THE TEMPERATURES OF SURFACES AND SURROtJNDINGS 

C BY THREE COLOR PYROMETRY. TWO OPTIONS ARE AVAILABLE: 

C 

C IFLAG = 1 IMPLIES THAT 

C THE SURFACE TEMPERATURE (TS) . THE EFFECTIVE SURROt^^TJING TEMPERATURE 

C (TF) , AND THE EMISSIVITY ARE KNOWti.' READING 1 (Rl), READING 2 (R2) , 

C AND READING 3 (R3) AND THE CORRESPONDING BRIGIITNESS 

C TEMPERATURES (BRTl) , (BRT2), AND (BRT3) ARE CALCUALTED. 

C 

C IFLAG = 2 IMPLIES THAT 

C Rl, R2. AND R3 (OR, BRTl, BRT2 , AND BRT3 ) ARE KNOWN. THE TS, TF, AND 

C THE EMISSIVITY ARE CALCULATED. 

C 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
•23456789012345678901234567890123456789012345678901234567890125456789012 

• REQUIRED INPUT 



IFLAG: IF IFLAG = 1 THEN CALCULATE Rl AND R2 . 

IF IFLAG = 2 THEN CALCUALTE TS AND TF. 
WLl = WAVELENGTH 1, MICRONS. 
WL2 = WAVELENGTH 2, MICRONS. 
WL3 = WAVELENGTH 1, MICRONS. 

EMIS = EMISSIVITY OF SURFACE AT WLl, WL2 , AND WL3 - 

TS = TEMPERATURE OF SURFACE, K. 

TF = EFFECTIVE SURROUNDING TEMPERATURE, K. 

Rl = READING 1 AT WAVELENGTH 1, W/M' 2 -MICRON . 

R2 = READING 2 AT WAVELENGTH 2, W/M*2 -MICRON . 

R3 = READING 3 AT WAVELENGTH 3, W/M'^ 2 • MICRON . 

BRTl = BRIGHTNESS TEMPERATURE FOR Rl, K. 

BRT2 = BRIGHTNESS TEMPERATURE FOR R2, K. 

BRT3 = BRIGHTNESS TEMPERATURE FOR R3 , K. 

ERR = ERROR RESULTING FROM OTHER SOURCES IN TERM OF BRIGHTNESS TEMPERATURE. K. 
EPS = CONVERGENCE CRITERION FOR NEWTON RAPHSON METHOD, W/M^2 . 
IMAX = MAXIMUM ALLOWED NUMBER OF ITERATIONS. 
DELTA = RELAXATION FACTOR. 



THE UNITS ON Cl AITO C2 ARE W-MICRONS*4 /M*2 AND MICRONS'K, RESPECTIVELY. 

2345678 90123 4 5678 90123 4 567 89012 345 678 9012 34567890 12 3 4 567 8 9012 3 4 56789012 
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PROGRAM TRICOL 

DIMEI'ISION ERR C 101 ) ,BRT1 (101) .BRT2 (101) , BRT3 (101) 
OPEN (14 , FILE= "rC8 .DAT* , STATUS= 'OLD' ) 
OPEN(18.rIbE='TC8.RES* ,STATUS=*NEW ) 

Cl=3 .742E8 

C2=l -439E4 
READ(14,') N, IFLAG, WL1,WL2,WL3 

WRITE(18,il) N, IFLAG, WLi,WL2. Wl,3 

11 FORMAT(lX,4X/5X, -N = ' , I2/5X. " IFI*AG = ■,I2/5X. 'WLl = •.F12.4, 

C ' MICRONS* /5X, •WL2 = •,F12.4.' MICROKS ' / 5X . * WL3 = ' ^12. 4, 

C • MICRONS ' ) 

IF (IFUAG.EQ-1) THEN 

READ (14,*) EMIS,TS,TF 

WF.ITEdS, 12) EMIS,TS,TF 
Wr.ITE( *, 12) EMIS,TS,TF 

12 FORMAT (IX, 4X/5X, -EMIS = * , F14 . 2 / 5X , ' TS = •.F14.1,' K' 
C /5X. 'TF = • ,F14 . 1. • K- ) 

Rl = FIAD(EMIS,TS,TF.WL1,C1,C2) 
R2 = RAD(EMIS,TS,TF.WL2,C1,C2) 
R3 = RAD(EMIS,TS,TF,WL3,C1,C2) 
raiT£(18,13) R1,R2,R3 
WRITE 13) R1,R2,R3 

13 FOPM.^T(lX.4X/5X, -Rl = •,£14.5.' W/M" 2 -MICRON ' /5X, ' R2 = 

C E14-S,' W/M-^S •MICRON ' /SX, 'R3 = ',E14.5,* W/M^^ 2 -MICRON • ) 

BRTll = BRT(R1,WL1.C1,C2) 
BRT12 = BRT(R2,WL2,C1,C2) 
BRT13 = BRT(R3,WL3,C1,C2) 
WT.ITE(18,17) BRT11,BRT12 , BRT13 
WRITE (-,17) BRTll, BRT12,BRT13 
17 FORMAT(:x,4X/5X, 'BRTl = •,F14.1,* K'/5X, •ERT2 = *. 

C F14.1.' KV5X,'BRT3= •,F14.1,* K' ) 

GO TO 155 
ELSE IF <IFUAG.E0.2) THEN 

REy^(14. •) EMISO,TSO,TFO, EPS, IMAX, DELTA 

WRITE (18,14) EMISO , TSO . TFO . EPS . IMAX , DELTA 
WRITE (*, 14) EMISO, TSO, TFO, EPS, IMAX, DELTA 

14 FORMAT ( IX, 4X/5X, 'EMISO = ■ , E14 . 2 / 5X, * TSO = •,F14.1,- K*/5X, 
C -TFO = •,F14.1,' K'/5X, -EPSs ',E14.5.' W/M^2' 

C /SX.'IMAX = I6/5X, 'DELTA = ',F12.4) 

DO 111 L=1,N 

READ ( 14 . • ) ERR(L) , BRTKL) , ERT2 (L) , BRT3 (L) 
111 CONTIinJE 

V^-RlTEtie. 1€) 
WRITER* . 1€) 

16 FORMAT (<X. -ERR' , 6X, ' BRTl ' , 6X. •BRT2' ,6X, •BRT3 ' , 6X, 'EMIS* , 6X, 

C 'TS* , 6X, *TF' , lOX, 'I-) 

DO 500 K=1,N 

D1IS=EMIS0 

TS^TSO 

TF=TFO 

Fl ^ EEL(BRT1 (K) ,WL1,C1,C2) 
R2 ' EBL(BRT2 (K) ,WL2,C1,C2) 
F3 = EEL(BRT3 (K) , WL3,C1,C2) 
El - P.1-RAD(EMIS,TS,TF, WL1,C1,C2) 
E: = P.2-RAD(EMIS,TS,TF,WL2,C1,C2) 
E3 = R3-RAD(EMIS,TS,TF. WL3,C1.C2) 
DO 20 1 = 1 , IMAX 

OEM = EMIS 
OTS = TS 

OTF = TF 

, All = EBL(TS.WLI,C1,C2)-EBL(TF,WL1.C1.C2) 
A12 = EMIS*DEDT(TS,WL1,C1,C2) 
A13 = (1-EMIS) •DEDT(TF, WL1,C1,C2) 
A21 = EBL(TS,V^T-2.C1,C2) -EBL(TF,WL2,C1.C2) 
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A2 2 = EMI S • DSDT (TS, , CI , C2 ) 
A23 = (i-EMIS) •DEDT(TF, WL2,C1,C2) 
A31 = EBL(TS,WL3,C1,C2) -EBL(TF,WL3,C1,C2) 
A32 = EMIS*DEDT(TS,WL3,C1,C2) 
A33 = (1-EMIS) *DEDT(TF, WL3 ,C1,C2) 
DETJl = A11-A22*A33+A12*A23-A31*A13*A21*A32 
DETJ2 = A13*A22*A31+A12*A21*A33'»-A11-A23-A32 
DETJ =DETJ1-DETJ2 
IF (DETJ.EQ.O.) THEN 
WRITE (18, 23) K.I 
WRITE(*,23) K,I 
23 FORMAT (IX, 4X, ' JACOB I AN IS SINGULAR AND ( K, I ) = 

C ( • ,14, • , • , 14, • ) • ) 

GO TO 100 

* 23 4 5 67890 12 34 567 89012 34 5678901 23 4567890123 45678 90 12 3 4 567 890 12 3 4 567 89012 
ENDIF 

DEI = B1*A22*A33+A12*A23*B3+A13*B2*A32 
DE2 = A13*A22*B3+A12*E2*A.33+B1*A23*A32 
DE = (DE1-DE2) /DETJ 

DSl = A11*B2*A33+B1*A23*A31+A13*A21-B3 
DS2 = A13*E2*A31+B1*A21*A33+A11»A23*B3 
OS = (DS1-DS2) /DETJ 

DFl = A11*A22*B3+A12»B2*A31+B1*A21*A32 
DF2 = B1*A22*A31+A12*A21*B3+A11*B2*A32 
OF = (DF1-DF2 ) /DETJ 
EMIS = OEM+DELTA*DE 
TS = OTS+DELTA*DS 
TF = OTF+DELTA*DF 

Bl ~ Rl-RAD(EMIS,TS,TF,vn.l,Cl,C2) 
B2 = R2-RAD(EMIS,TS,TF,WL2,C1,C2) 
B3 = R3-RAD(EMIS,TS,TF>7L3,C1,C2) 
TEST = SQRT(B1*B1+B2*B2+B3*B3) 
IF (TEST.LE.EPS) GO TO 99 

20 CONTINUE 
ELSE 

WRITE (•,21) 
WRITEdS. 21) 

21 FORMAT (IX, 4X, • INAPPROPRIATE VALUE OF IFLAG. ' ) 
GO TO 155 

ENDIF 

WRITE{18,15) K,I 
WRITE t VIS) K,I 
15 FORMAT ( IX, 4X, 'IMAX EXCEEDED WITHOUT CONVERGENCE AND 
C (K,I) = (M4, * , ' ,14, •) •) 
GO TO 100 
99 CONTINUE 

WRITE(*,19) ERR(K) ,BRT1(K) ,BRT2(K) , BRT3 <K) .EMIS,TS,TF. I 
WRITEdS, 19) ERR(K) ,BRT1(K) ,BRT2(K) , BRT3 (K) ,EMIS,TS,TF. I 
19 FORMAT ( IX , F8 . 1 , 2 X , F8 . 1 , 2X , F8 . 1 , 2X, F8 . 1 , 2X , F6 . 4 , 2X. F8 . 1 , 
C 2X.F8.1.2X,I6) 
100 CONTINUE 
500 CONTINUE 
155 CONTINUE 
CLOSE (18) 
PAUSE 
STOP 
END 

* 2 345678901234 56789012 34 567 8901234 5678901 234 567890123 4 56789012 3 4 5678 9012 



FUNCTION rad(e:mis,ts,tf,wl, C1,C2 ) 

RAD = (1 . -EMIS) •EBL(TF,WL,C1 , C2) +EMIS*EBL(TS,WL, CI ,C2 ) 

RETURN 

END 
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FUNCTION EBL(T,WL,C1.C2) 

EBL = C1/(WL-*5*(E:XP{C2/(WL»T))-1.)) 

RETURN 

END 



FUNCTION BRT(E,WL.C1,C2) 

BRT = C2/ (WL*ALOG( {CI/ (R*WL**5) +1 - ) ) ) 

RETURN 

END 



FUNCTION DEDT(T, WL,C1,C2) 

DEDT = Cl*C2*E:xr(C2/ (WL*T) ) / (WL**6*T-T* (EXP(C2/ (WL*T) ) -1 . ) **2) 

RETURN 

END 



*23456789012345678901234S67890123456789012 34567890123456789012 34 56789012 



(2) A sample input file: 

11 
1 

0. 8D0 

1 . 6D0 
2.2dO 

0. 4DO 
1300. DO 
1500. DO 

1. D-1 
1000 
1,D0 



-5.00 


1450 . 


,3 


1434 


.1 


1424 


.0 


-4.00 


1450. 


,3 


1434 


.1 


1425 


.0 


-3.00 


1450. 


3 


1434 


.1 


1426 


.0 


-2.00 


1450. 


3 


1434 


.1 


1427 


.0 


-1,00 


1450. 


3 


1434 


.1 


1428 


.0 


0.00 


1450. 


3 


1434 


.1 


1429 


.0 


1.00 


1450. 


3 


1434 


.1 


1430 


.0 


2.00 


1450. 


3 


1434 


.1 


1431 


.0 


3.00 


1450. 


3 


1434 


.1 


1432 


.0 


4.00 


1450. 


3 


1434 


.1 


1433 


.0 


5.00 


1450. 


3 


1434 


.1 


1434. 


.0 



(3) A sample output file: 

N = 11 

IFIAG = 2 

WLl = .8000 MICRONS 

WL2 = 1.6000 MICRONS 

WL3 = 2.2000 MICRONS 

EMISO ~ 0.40E+00 

TSO = 1200.0 K 

TFO = 1600.0 K 

EPS = 0.1000OE*00 W/M*2 

IMAX = 1000 
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DELTA ^ 1.0000 



ERR 




BRTl 




3RT2 




BRT3 




-5 


.0 


1450 


. 3 


1434 


. 1 


1424 


. 0 


-4, 


.0 


1450 


.3 


1434 


. 1 


1425 


. 0 


-3 . 


. 0 


1450 


. 3 


1434 


. 1 


1426 


- 0 


-2 . 


, 0 


1450 


. 3 


1434 


. 1 


1427 , 


. 0 


-1. 


,0 


1450. 


.3 


1434 . 


. 1 


1428, 


,0 




0 


1450 


.3 


1434 . 


. 1 


1429, 


.0 


1. 


0 


1450. 


.3 


1434 . 


, 1 


1430. 


,0 


2 . 


0 


1450. 


3 


1434 . 


1 


1431 . 


0 


3 . 


0 


1450. 


3 


1434- 


1 


1432. 


0 


4 . 


0 


1450. 


3 


1434 . 


1 


1433 - 


0 



JACOBIAN IS SINGULAR AND ( K, I ) = 
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EMIS 


TS 




TF 




I 


1360 


885. 


. 3 


1467 


. 6 


6 


1445 


979. 


.0 


1468 


.7 


6 


1602 


1060. 


. 8 


1470, 


.8 


6 


1901 


1139. 


. 5 


1474 . 


. 6 


5 


2525 


1219 . 


3 


1482 . 


, 2 


5 


3994 


1299 . 


0 


1499. 


9 


5 


6979 


1369 . 


c 


1548 . 


6 


4 


9280 


1408 . 


3 


1666- 


0 


11 


9879 


1424 . 


1 


1887, 


9 


13 


9987 


1430. 


9 


2374 . 


0 


10 


( 11, 


7) 
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CLAIMS 

What is claimed is : 

1 . A method of measuring the temperature of a surface, said method 
comprising: 

5 providing a surface; 

selecting a first wave length; 
providing a first pyrometer; 

adapting said first pyrometer to sense radiation of said first wave length; 

selecting a second wave length; 
10 providing a second pyrometer; 

adapting said second pyrometer to sense radiation of said second wave length; 

focussing said first pyrometer on said surface and obtaining a first reading; 

focussing said second pyrometer on said surface and obtaining a second reading; 

determining an emissivity of said surface; and 
15 determining said temperature by utilizing said first reading; said second reading and 
said emissivity. 

2. The method of claim 1 further comprising positioning said surface within 
an enclosure prior to said focussing of said first pyrometer on said surface. 

20 

3. The method of claim 1 further comprising selecting said first wave length 
to be within the range of 0.5 to 5 microns. 

4. The method of claim 1 further comprising selecting said second wave 
25 length to be within the range of 0.5 to 5 microns. 

5. The method of claim 1 further comprising determining said emissivity 

by: 

insulating said surface; 
30 heating said surface; 

focussing a third pyrometer on said surface; 

obtaining readings from said third pyrometer as said surface cools; and 
determining said emissivity by utilizing said readings from said third pyrometer. 
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6. The method of claim 1 further comprising adapting said first pyrometer 
to sense radiation of said first wave length by using a filter. 

7. The method of claim 1 further comprising adapting said second 
5 pyrometer to sense radiation of said second wave length by using a filter. 

8. The method of claim 1 wherein said first pyrometer and said second 
pyrometer are one and the same pyrometer. 

10 9. The method of claim 1 wherein said first wave length is selected to 

minimize the effects of particulate material in the line of sight of said first 
pyrometer to said surface, 

10. The method of claim 1 wherein said second wave length is selected to 
15 minimize the effects of particulate material in the line of sight of said second 

pyrometer to said surface. 

11. The method of claim 1 further comprising: 
selecting a third wave length; 

20 providing a third pyrometer; 

adapting said third pyrometer to sense radiation having said third wave length; and 
focussing said third pyrometer on said surface and obtaining readings. 
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